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ABSTRACT 


The  reflection  of  a  spherical  acoustic  wave  by  a  flat  Impedance 
plane  is  a  deceptively  simple  problem.  There  is,  in  fact,  no  exact 
closed-form  solution.  In  practice,  approximate  methods  are  employed 
to  predict  the  received  pressure,  the  most  popular  of  which  relies  on 
the  use  of  the  classical  plane  wave  reflection  coefficient,  which 
proves  to  be  deficient  for  grazing  angles  of  incidence. 

The  primary  objective  of  this  study  was  to  derive  a  more  exact 

a.c0usftc,  , 


solution  to  the  problem  of^ point-source  propagation  over  a  locally- 
reacting^0  Impedance  or ^’'grounf'*5 plane.  This  objective  was  met  with 
the  derivation  of  an  asymptotic  series  solution.  One  of  the  most 
important  features  of  this  solution  is  that  higher-order  terms  can  be 
calculated  from  preceding  terms  in  the  series  by  the  use  of  recursion 
formulae,  also  derived  here.  Comparing  data  predicted  from  this 
solution  with  that  from  a  numerical  integration  of  the  exact 
expression  showed  the  asymptotic' series  to  be  extremely  accurate,  even 
for  very  low  values  of  the  parameter  kR.  As  expected,  the  plane  wave 
solution  often  showed  major  deviations  from  the  exact  integral 
solution.  ~~~? 


A  secondary  goal  was  to  incorporate  the  new  propagation  solution 
into  a  barrier  model  so  that  ground  reflections  in  addition  to  edge 
diffraction  could  be  accounted  for.  Only  the  first  term  in  the 
asymptotic  ground  propagstlon  solution  was  used  for  this  purpose,  hs 
it  was  shown  to  be  sufficiently  accurate  for  many  practical  cases. 
Thus,  an *"£dge-Plus-Image8* ^ barrier  diffraction  model  was  developed  in 


-^-the  second  phase  of  this  study.  The  results  of  preliminary 

sensitivity  tests  reported  here  are  very  encouraging,  and  indicate 
that  the  barrier  model  should  afford  a  higher  degree  of  accuracy  than 
available  with  similar  models  employing  the  plane  wave  reflection 
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CHAPTER  I 


INTRODUCTION 

The  problem  of  spherical  wave  acoustic  propagation  over  an' 
Impedance  plane  continues  to  attract  attention  within  the  scientific 
community,  although  the  results  of  the  first  major  theoretical  studies 
on  the  subject  appeared  in  the  literature  over  thirty  years  ago.  Much 
consideration  is  also  being  given  to  the  effects  on  the  sound  field 
when  a  diffracting  barrier  is  located  between  the  source  and  receiver. 
At  first  glance,  such  strong  research  Interest — a  reliable  indication 
that  these  problems  are  far  from  "solved" — may  be  surprising, 
especially  in  view  of  the  fact  that  many  of  the  underlying 
mathematical  concepts  were  well-developed  over  a  century  ago  in  the 
field  of  optics,  and  further  refined  in  the  early  1900" s  in  the  field 
of  electromagnetics.  However,  in  response  to  more  contemporary 
research  needs — particularly  those  of  underwater  acoustics  as  well  as 
community  and  aircraft  noise  control — a  re-evaluation  of  previous 
results  has  become  necessary. 


1.1  Background  and  Statement  of  the  Problem 


Figure  1.1  (a)  and  (b)  shows  the  geometry  of  the  overall  problem, 
where  sound  propagates  from  the  source  S  to  the  receiver  R  over  a 
diffracting  plane  B  on,  and  perpendicular  to,  an  Impedance  plane  G. 
Although  delineated  for  clarity,  the  Impedance  plane  is  actually 
infinite  in  both  of  its  dimensions,  while  the  diffracting  plane  is 
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Infinite  in  on*  dimension  and  finite  in  the  other.  The  Infinite 
entente  ere  represented  by  snail  arrows  in  Figure  1.1  (a).  Without 
loes  of  generality,  the  inpedance  plane  G  can  be  assumed  to  be 
"horizontal”  and  therefore  referred  to  as  the  "ground  plane,"  while 
the  diffracting  plane  B  can  be  assumed  "vertical,"  of  finite  "height," 
infinitely  "long,"  and  referred  to  ae  the  "barrier."  The  term 
"barrier"  is  appropriate  for  the  diffracting  plane  since  it  will 
generally  be  assumed  in  this  study,  as  is  usually  the  case  in 
practice,  that  the  receiver  is  in  its  acoustic  "shadow"  relative  to 
the  source. 

The  Impedance  of  the  ground  plane  may  take  different  values  on 
the  source  and  receiver  aides  of  the  barrier  (Z^  and  in  Figure  1.1 
[b]);  however,  it  will  be  assumed  constant  on  either  side.  A 
fundamental  constraint  here  is  that  the  impedance  of  the  ground  plane 

does  not  depend  on  the  angle  of  the  incident  acoustic  energy,  and  so 

✓ 

it  may  be  fully  described  by  specifying  only  its  value  at  normal 
incidence.*' 

The  flow  of  acoustic  energy  from  the  source  to  the  receiver  may 
be  regarded  as  proceeding  by  way  of  four  distinct  paths,  shown  as  rays 
with  exaggerated  spacing  in  Figure  1.1  (a).  Thus,  there  is  (i)  a 
direct  path  over  the  barrier  without  ground  interaction  (path  S8R  in 
Figure  1.1  [b]);  (ii)  a  reflect ed-dlrect  path  (SP^BR)  for  energy  that 

1  This  is  the  same  ae  saying  that  the  surface  will  not  support 
acoustic  particle  velocities  at  angles  other  then  normal,  or  that  the 
surface  is  "locally  reacting."  It  is  becoming  common  to  use 
"impedance  plane"  to  imply  these  characteristics. 


la  first  reflected  from  the  ground  on  the  source  side  of  the  barrier 
and  then  diffracted  over  the  barrier  to  the  receiver;  (111)  a  direct- 
reflected  path  (SBP2R),  where  the  ground  reflection  is  on  the  receiver 
side  of  the  barrier;  and  (lv)  a  reflected-reflected  path  (SP^BP^)  for 
energy  first  reflected  from  the  source-side  ground,  then  diffracted  by 
the  barrier,  and  finally  reflected  to  the  receiver  from  the  receiver- 
side  ground. 

The  overall  problem  of  a  barrier  on  a  ground  plane  splits  quite 
naturally  Into  two  parts:  propagation  over  an  impedance  plane  and 
diffraction  by  a  barrier.  The  former  is  illustrated  In  Figure  1.2  (a) 
and  (b).  This  deceptively  simple  problem  has  no  closed-form  solution. 
Complexities  derive  from  the  fact  that  the  spherical  waves  from  the 
point  source  are  not  easily  matched  to  the  planar  boundary  conditions 
dictated  by  the  Impedance  plane. 

It  is  with  regard  to  this  problem  that  the  first  major 
deficiencies  in  classical  theory  became  apparent.  Older  studies 
assumed  that  the  incident  acoustic  energy  was  in  the  form  of  plane 
waves,  for  which  the  problem  Is  amenable  to  solution.  However,  since 
plane  waves  are  a  fiction  of  mathematics  (implying  a  source  of  energy 
"infinitely  far"  away),  the  resulting  solution  becomes  Inaccurate  when 
applied  to  more  realistic  source-receiver  configurations.  Moreover, 
for  angles  of  Incidence  that  approach  grazing,  the  plane  wave  theory 
Is  not  useful  at  all.  Bence,  modern  researchers,  working  on  problems 
where  grazing  incidence  Is  the  rule,  have  been  forced  to  take  into 
account  the  sphericity  of  the  incident  wavefronts. 
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Figure  1.2  The  Geometry  of  the  Ground  Propagation  Problem 
(e)  perspective  view,  (b)  aide  view. 


The  second  part  of  the  problem  under  study,  the  diffraction  by 
the  barrier,  is  illustrated  in  Figure  1.3  (a)  and  (b).  The  case  for 
which  the  barrier  extends  to  infinity  in  the  two  directions  indicated 
by  the  small  arrows  (thus,  a  "half  plane"  or  “semi-infinite"  plane) 
has  received  considerable  attention  in  the  past.  Yet  it  is  clear  that 
a  barrier  located  on  a  ground  plane  as  in  Figure  1.1  does  not  fit  such 
a  case,  and  to  use  classical  results  without  modification  to  describe 
the  diffraction  process  here  would  be  imprecise.  The  effects  of  the 
ground  plane  must  be  taken  into  account,  and  these  effects  become  more 
pronounced  as  the  barrier  height  decreases. 

The  research  to  be  described  in  this  thesis  has  retained,  and 
expanded  upon,  many  of  the  constructs  of  classical  half-plane 
diffraction  theory  in  formulating  a  model  for  the  overall  problem 
depicted  in  Figure  1.1;  however,  a  scheme  has  been  devised  for  fusing 
these  constructs  with  the  results  of  an  original,  more  exact  solution 
to  the  ground  propagation  problem  of  Figure  1.2.  Predictions,  based 
on  this  approach,  to  acoustic  propagation  problems  involving  a 
diffracting  barrier  located  on  a  ground  plane  should  result  in  a 
higher  degree  of  accuracy  than  available  from  current  models. 

1.2  Review  of  Previous  Investigations 
The  composite  problem  of  acoustic  propagation  over  a  ground  plane 


in  the  presence  of  a  diffracting  barrier  has  been  addressed  only 
recently.  But,  as  mentioned  in  the  preceding  section,  the  components 
parts  of  the  problem—nauely,  reflection  at  a  plane  Interface  and 


(b) 


Figure  1.3  The  Geonetty  of  eh*  Half-Plane  Diffraction 

Problem.  (a)  perspective  view,  (b)  side  view 


diffraction  by  a  half  plane — are  classical  and  have  been  well-studied 
separately. 


1.2.1  The  Problem  of  Propagation  Over  an  Impedance  Plane 

The  origins  of  the  formal  reflection  problem  were  in  the  field  of 

optics.  Thomas  Young  (1773-1829)  and  Augustin  Fresnel  (1788-1827) 

performed  the  earliest  work  on  plane  waves  at  perpendicular  incidence 

to  an  interface.  In  the  mld-1800's,  Simeon  Poisson,  George  Green, 

Lord  Rayleigh,  and  Herman  Helmholtz  researched  the  analogous 

2 

acoustical  reflection  problem,  including  oblique  Incidence. 

Sommerfeld  (1909,  1926),  working  in  the  field  of  electro¬ 
magnetics,  was  the  first  to  solve  the  spherical-wave  reflection 
problem,  stated  as  a  dipole  source  on  a  finitely-conducting  earth. 

Weyl  (1919)  re-formulated  the  problem  by  modelling  the  radiation  from 
a  point  source  located  above  the  earth  as  a  superposition  of  an 
infinite  number  of  elementary  plane  waves,  propagating  in  different 
(complex)  directions.  In  the  derivation,  each  component  wave  is 

3 

reflected  according  to  the  classical  plane  wave  laws. 


A  good  discussion  of  such  historical  foundations  may  be  found 
in  Lindsay  (1972). 

3 

Although  the  final  forms  of  the  Sommerfeld  and  Weyl  solutions 
are  very  similar,  the  former  contains,  while  the  latter  lacks,  an 
explicit  "surface  wave”  term.  Much  controversy  followed  on  this 
subject  (some  of  which  continues  today  for  the  acoustical  problem), 
the  nature  of  which  can  be  gleaned  from  Norton  (1937;  pp.  1193-1195, 
pp.  1234-1236)  or  Stratton  (1941,  pp.  573-587). 
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The  Sommerfeld  and  Weyl  integral  formulations  were  not  conducive 
to  numerical  calculations  or  physical  insight  until  Van  der  Pol  (1935) 
applied  several  ingenious  substitutions  that  simplified  certain 
Integrals  appearing  in  the  derivations.  Norton  (1936,  1937)  expanded 
upon  these  and  other  results  from  Van  der  Pol  and  Niessen  (1930)  and, 
with  the  aid  of  equations  by  Vise  (1931),  generated  the  most  useful 
results  up  to  that  time.  His  final  expressions  for  the  field,  though 
approximate,  were  in  terms  of  the  complex  complementary  error 
function,  for  which  numerical  values  could  be  calculated.  Although 
scattered  research  continued  on  the  electromagnetic  reflection 
problem,  it  is  fair  to  say  that  the  inquiry  was  essentially  completed 
with  Norton's  publications. 

The  acoustical  problem  of  spherical  wave  refection  was  first 
attacked  by  Rudnlck  (1947),  who  relied  heavily  on  the  electromagnetic 
theories  of  Van  der  Pol  and  Norton.  He  adapted  and  applied  these 
results  to  the  problem  of  determining  the  sound  field  from  a  point 
source  in  the  vicinity  of  a  plane  Interface  between  two  fluid  media. 

In  subsequent  papers  with  Lawhead  (Lawhead  and  Rudnlck,  1951a,  1951b), 
Rudnlck  extended  his  solution  to  apply  to  a  point  source  located  above 
an  Impedance  plane.  Approximations  analogous  to  those  made  by  Norton 
also  enabled  Rudnlck  and  Lawhead  to  reduce  the  integrals  appearing  in 
the  Van  der  Pol  formulation  to  the  error  function  form,  from  which 
they  obtained  and  plotted  numerical  results.  About  the  same  time, 
Ingard  (1951)  applied  tfeyl's  method  to  the  impedance  plane  problem  and 
obtained  results  very  slnllar  to  Rudnlck's,  although  the  methods  of 
analysis  ware  quite  different. 
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More  recent  theoretical  studies  have  expanded  upon  the 

foundations  laid  by  Rudnick,  Lauhead,  and  Ingard.  Paul  (1957)  studied 

the  two-media  problem,  applying  the  method  of  steepest  descent  as 

modified  by  Bands  and  Wesley  (1953),  to  the  Sommerfeld-Rudnick  theory 

and  obtained  two  asymptotic  solutions,  one  valid  near  the  interface 

and  the  other  valid  in  the  region  above  the  source.  Wenzel  (1974) 

analyzed  the  impedance  plane  problem  via  Green's  functions  and  contour 

Integration,  and  obtained  Integrals  to  which  he  applied  asymptotic 

techniques  for  various  extremes  of  the  surface  impedance.  The 

"surface  wave”  was  attributed  much  importance  by  Wenzel  and  an 

explicit  term  appeared  in  his  solution.  Chien  and  Soroka  (1975,  1980) 

applied  Bands'  (1966)  double  saddle  point  method  of  integration — the 

"method  of  subtraction  of  the  pole" — to  derive  an  asymptotic 

expression  for  the  total  field,  which  contained  the  complementary 

2 

error  function.  Their  solution  is  explicit  through  terms  in  1/R  , 

where  R  is  the  distance  from  the  source  image  in  the  plane  to  the 
receiver.  Thomasson  (1976,  1977)  employed  the  method  of  steepest 
descent,  along  with  several  variable  transformations,  and  also 
expressed  his  solution  in  terms  of  the  error  function.  His  result  is 
very  similar  to  the  first  term  of  Chien  and  Soroka's  full  solution. 

In  a  subsequent  effort  (concurrent  with  the  present  study),  Thomasson 
(1980)  generated  an  asymptotic  expansion  of  his  exact  integral 
representation.  He  also  gave  recursion  relations  for  the  coefficients 
in  this  expansion  to  enable  higher-order  terms  to  be  derived  or,  more 
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readily i  calculated  on  a  digital  computer. 

In  addition  to  the  principal  studies  above,  there  have  been  many 
other  investigations — some  theoretical,  but  most  experimental — on  the 
subject  of  sound  propagation  above  an  impedance  plane  over  the  past 
decade.  The  following  selected  listing  may  serve  as  a  guide  for 
further  information:  Attenborough  (1982);  Attenborough,  Hayek,  and 
Lawther  (1980);  Briquet  and  Filippi  (1977);  Butov  (1981);  Chessell 
(1977);  Delany  and  Bazley  (1970);  Donato  (1976a,  1976b);  Embleton, 
Piercy,  and  Olson  (1976);  Filippi  and  Habault  (1978);  Habault  (1980); 
Habault  and  Filippi  (1981);  Hayek,  Lawther,  Kendig  and  Simowitz 
(1978);  Hayek,  Attenborough,  and  Lawther  (1980);  Hayek,  Lawther,  and 
Tate  (1980);  Lawther,  Hayek,  Tate,  and  Nobile  (1980);  Naghieh  and 
Hayek  (1981);  Pao  and  Evans  (1971);  Rassmussen  (1982);  and  Van  faoorhem 
(1975). 


1.2.2  The  Problem  of  Diffraction  by  a  Barrier 


With  its  origins  in  optics  and  its  development  in 
electromagnetics,  the  study  of  wave  diffraction  around  obstacles  has 
experienced  a  new  surge  of  interest  in  the  field  of  acoustics, 
primarily  among  scientist  working  in  community  noise  control  and 
underwater  acoustics.  The  body  of  literature  on  the  subject  of 


diffraction  is  immense,  even  when  limited  to  studies  of  the  half 


Although  not  realized  by  Thomasson,  it  has  been  shown  in  the 
course  of  the. present  study  that  Chien  and  Soroka's  full  solution  (to 
terms  in  1/RZ'  could  be  extracted  Identically  from  the  first  few  terms 
in  Thomasson's  asymptotic  solution. 


plane.  Thus,  only  a  general  survey  will  be  given  here. 

Thomas  Young,  In  1800,  was  the  first  to  critically  study  the 
diffraction  problem,  but  his  assertions  that  the  phenomenon  was  local 
to  the  diffracting  edge  were  rejected  by  Fresnel  in  1815.  Instead, 
Fresnel  employed  Huygens'  principle  and  the  wave  nature  of  light  to 
develop  a  mathematical  theory  which  became  the  foundation  for  more 
modern  research  in  diffraction.  Helmholtz,  in  1859,  and, 
independently,  Kirchoff,  in  1882,  derived  their  well-known  integral 
equation  relating  the  field  at  an  arbitrary  point  in  a  region  to  the 
value  of  the  field  and  its  derivatives  on  the  bounding,  or 
diffracting,  surface. 

Modern  diffraction  theory  has  its  roots  in  the  research  of 
Sommerfeld  (1896),  whose  rigorous  approach  using  an  extension  of  image 
theory  yielded  an  exact  integral  expression  for  the  diffraction  of  a 
plane  wave  by  a  rigid  (or,  alternatively,  perfectly  soft)  naif  plane. 
Carslaw  (1898)  continued  Sommerfeld's  research  on  the  half  pi  -  ~  and 
derived  a  solution  for  a  point  source.  MacDonald  (1901,  911) 
addressed  the  problem  of  wedge  diffraction — for  both  a  point  and  a 
line  source — using  the  classical  "separation  of  variables"  technique 
and  expressed  his  solution  as  an  infinite  series  of  appropriate 
eigenfunctions.  He  showed  how  this  series  could  be  re-written  in 
integral  form  and  noted,  for  the  particular  case  when  the  wedge 
"collapsed"  to  a  straight  edge,  that  this  integral  solution  was 
identical  to  that  given  previously  by  Carslaw.  MacDonald  went  on  to 
derive  an  asymptotic  form  for  his  exact  solution  in  terms  of  Fresnel 
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Integrals,  a  form  which  is  utilized  often  in  modern  diffraction 
studies. 

Redfearn  (1940),  recognizing  "the  practical  Importance  of  the 
problem  of  reduction  of  noise  at  a  point  due  to  an  unavoidable  source 
of  sound  at  a  second  point,”  (p.  273)  was  the  first  to  seek  a 
simplified,  engineering  approach  and  solution  to  the  edge  diffraction 
problem  for  the  acoustical  case.  Drawing  analogies  from  heat 
conduction  theory,^  he  put  Sommerf eld's  solution  into  a  form  for  a 
point  source,  discarded  the  second  of  two  integral  terms  (which,  in 
fact,  represents  the  source  image  in  the  diffracting  plane),  assumed 
the  gradient  of  the  field  is  zero  on  the  faces  of  the  barrier  (an 
assumption  first  proposed  by  Kirchoff),  and  finally,  generated  a 
solution  in  terms  of  a  single  geometrical  parameter.  Redfearn' s  paper 
seems  to  have  attracted  little  attention  at  the  time  of  its 
publication,  but  it  was  re-discovered  in  the  early  1960's  when  it 
served  as  the  impetus  for  more  modern  single-number  prediction  schemes 
for  attenuation  of  community  noise  by  barriers. 

In  response  to  the  advances  in  radio  communications  in  the  late 
1940's  and  early  1950's,  researchers  began  looking  for  practical  means 
of  obtaining  accurate  numerical  results  for  the  electromagnetic 
diffraction  problem.  Thus,  the  Weiner-Hopf  technique  became  popular 
(Heins,  1956;  Rawlins,  1975;  Senior,  1952)  and  so-called  variational 
techniques  were  applied  to  the  problem  (Levine  &  Schwinger,  1948, 

^  Redfearn  makes  reference  to  "Carslaw,  'Conduction  of  Heat',” 
but  gives  no  citation. 


1949).  The  duel-integral  method,  which  yields  two  integral  equations 
valid  in  separate  regions  of  the  complex  plane,  was  also  employed 
(Clemmow,  1966;  Jones,  1952).  Although  applied  with  success,  none  of 
these  methods  was  especially  convenient. 

The  breakthrough  that  facilitated  practical  solutions  to  many 
common  diffraction  problems  was  due  to  Keller  (1958,  1962)  and  his 
work  on  the  Geometrical  Theory  of  Diffraction  (GTD) .  Essentially,  the 
GTD  method  extends  classical  geometric  acoustics  to  Include  not  only 
reflected  sound  rays  but  diffracted  rays  as  well.  The  mathematical 
description  of  the  diffracted  rays  is  in  terms  of  a  "diffraction 
coefficient"  deduced  from  an  approximation  to  a  "canonical" 
diffraction  solution.  For  example,  in  the  treatment  of  a  half  plane 
Keller  used  the  first  term  in  the  asymptotic  expansion  of  Sommerfeld's 
exact  plane  wave  solution  to  derive  the  diffraction  coefficient.  His 
resulting  GTD  solution,  however,  is  not  valid  in  the  transition 
regions  of  the  geometric  shadow  boundaries.  Recognizing  this 
shortcoming,  Kouyoumjian  and  Fathak  (1974)  incorporated  previous 
asymptotic  results  from  Pauli  (1938)  and  Oberhettlnger  (1956)  into 
Keller's  Geometrical  Theory  of  Diffraction  to  generate  a  uniform 
asymptotic  solution,  in  terms  of  Fresnel  integrals,  which  remained 
valid  in  the  transition  regions  of  the  shadow  boundaries. 

Concerning  the  more  complex  problem  of  the  half  plane  with 
arbitrary  surface  impedances,  Malyuzhlnets  (1955,  1962)  and 
Malyuzhinets  &  Tuzhlln  (1970)  employed  generalized  Fourier  transforms 
to  obtain  the  first  rigorous  solution.  However,  this  was  left  in 
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Integral  form,  and  special,  cumbersome  functions  were  needed  for  its 
evaluation.  Recently,  Kendig  (1977)  derived  a  more  useful,  closed- 
fora  solution  for  this  problea  (as  well  as  for  the  cases  of  rigid  and 
soft  half-plane  faces — see  Kendig  &  Hayek,  1981)  using  function 
theoretic  aethods.  Kendig's  results  were  expressed  in  teras  of 
Fresnel  integrals,  which  are  relatively  siaple  to  evaluate 
nuaerically.  A  comprehensive  research  effort  by  the  Applied  Research 
Laboratory  of  The  Pennsylvania  State  University  (Hayek,  et  al.,  1978) 
consolidated  Keller's  Geoaetrlcal  Theory  of  Diffraction,  Kouyouajian's 
and  Pathak's  uniform  asymptotic  results,  and  Kendig's  diffraction 
solution,  to  derive  very  accurate,  manageable  expressions  for  the 
barrier  problea.  These  expressions  serve  as  the  basis  for  the 
diffraction  work  of  the  present  study  (see  also  Hayek,  1982;  Lawther  & 
Hayek,  1978;  Lawther  et  al.,  1980;  and  Hayek  &  Nobile,  1981). 

Excellent  tutorial  and  review  material  on  the  subject  of 
diffraction,  including  extensive  bibliographies  may  be  found  in 
Pierce  (1974);  Skudrzyk  (1971);  Bouwkaap  (1954);  and  Bowman,  Senior,  & 
Uslenghl  (1969). 


1.2.3  The  Combined  Problea  of  Diffraction  by  a  Barrier  on  a  Ground 


Plane 


Due  to  the  complicated  boundary  geometry,  there  is  no  exact, 
closed-fora  solution  to  the  problea  of  point  source  diffraction  by  a 
thin,  planar  barrier  located  on  a  flat,  impedance-covered  ground 
plane.  Neither  is  there  an  exact  Integral  fora  of  the  solution  that 


can  be  evaluated  numerically.  Even  when  both  the  barrier  and  the 
ground  plane  are  perfectly  rigid,  an  exact  solution  can  only  be 
formulated  (Nobile  &  Hayek,  1981)  in  terms  of  an  Integral  over  an 
Infinite  summation  of  Mathleu  functions,  and  to  extract  any  numerical 
results  from  this  can  be  most  difficult.  Lacking,  then,  a  general 
solution  to  which  mathematical  approximations  could  be  applied,  the 
method  of  attack  In  most  studies  has  been  one  of  making  physical 
assumptions  about  the  sound  field.  The  fact  that  these  assumptions 
differ  from  one  author  to  another  is  principally  responsible  for  the 
variety  of  solutions  offered  in  the  current  literature.  Sometimes 
very  broad  and  admittedly  crude  assumptions  produce  theoretical 
predictions  that  compare  favorably  with  experiments.  However,  an 
Increase  in  accuracy  is  almost  always  attainable  using  fewer,  or  less 
broad,  assumptions.  The  price  for  this  la  much  more  complexity  In  the 
mathematics  and  often  more  difficulty  in  extracting  numerical  results. 
Although  there  have  been  many  experimental  studies  on  this  topic,  the 
theoretical  work  Is  sparse. 

Redfearn  (1940),  mentioned  above,  spoke  in  terms  of  applying  his 
results  to  actual  barriers  erected  on  the  ground,  but  he  proposed  no 
method  for  handling  the  ground  reflections.  Fehr  (1951),  on  the  other 
hand,  specifically  addressed  the  ground  probleo,  but  assumed  that  the 
ground  was  perfectly  rigid  and  that  both  source  and  receiver  lay 
directly  on  it;  consequently,  his  approach  simply  became  a  matter  of 
accounting  for  the  pressure  doubling  by  the  ground. 


Maekawa  (1965,  1968)  vas  tha  flrat  to  conduct  a  comprehensive  set 
of  experiments  with  the  expressed  purpose  of  obtaining  a  large  data 
base  for  comparison  with  existing  theories.  He  found  that  his 
experimental  points  clustered  rather  close  to  Redfearn's  predicted 
values,  and  by  fitting  a  line  through  the  points  plotted  against  the 
so-called  Fresnel  number,  Maekawa  offered  a  general  design  curve  for 
calculating  the  attenuation  expected  from  a  half-plane  barrier. 

Maekawa  went  on  to  propose  a  simple  scheme  based  on  his  empirical 
curve  to  account  for  the  effects  of  (rigid)  ground  reflections  on  the 
receiver  side  of  the  barrier  (he  assumed  the  source  was  still  directly 
on  the  ground).  Rathe  (1969)  re-phrased  Maekawa's  results  into  a  form 
that  easily  yields  engineering  predictions.  Kurze  and  Anderson  (1971) 
showed  that  the  prediction  schemes  of  Maekawa  and  Rathe  also  remain 
valid  for  oblique  incidence  and  additionally  proposed  an  equation  that 
closely  describes  the  empirical  curve. 

The  first  theoretical  studies  that  addressed  the  problem  of  a 
finite-impedance  ground  plane  and  that  took  into  account  ground 
reflections  from  the  source  side  of  the  barrier  in  addition  to  those 
from  the  receiver  side  were  done  by  Lindblad  (1970)  and  Jonasson 
(1972).  Lindblad  set  up  a  numerical  integration  of  elementary  sources 
(Helmholtz-Huygens  Integral  approach^)  over  the  infinite  plane  above 
the  barrier.  The  major  assumptions  in  his  approach  are  that  all 
ground  reflections  behave  in  accordance  with  plane  wave  theory,  and 
that  the  fields  above  the  barrier  are  precisely  what  they  would  be  In 

6  For  example,  see  Skudrzyk  (1971,  chap.  23). 


18 


the  absence  of  the  diffracting  barrier  (the  so-called  "Kirchoff 
assumption") .  Jonasson  assuaed  that  the  ground  was  locally-reacting 
and  that  the  barrier  was  totally  absorptive.  The  latter  condition 
enabled  hla  to  Ignore  the  second  tern  in  the  MacDonald  solution,  which 
Jonasson  employed  to  describe  the  diffraction  effects.  The  problem 
was  treated  as  a  simple  superposition  of  four  separate  half-plane 
diffraction  problems  (see  Figure  1.4),  where  the  "strengths”  of  the 
ground-reflected  rays  were  derived  from  spherical-weve  propagation 
theory  (Ingard,  1951). 

Thonasson  (1978)  correctly  pointed  out  that  neither  Llndblad's 
nor  Jonasson' a  solution  is  valid  for  small  barrier  heights.  In  the 
former,  the  plane-wave  reflection  theory  falls  for  the  near-grazing 
rays  required  to  construct  the  fields  above  the  barrier.  In  the 
latter,  the  Inherent  assumption  that  each  ray  "sees"  a  semi-infinite 
barrier  becomes  unsound.  Thonasson  offered  an  Integral  solution  that 
is  more  accurate  for  short  (and  also,  for  the  first  time,  finite- 
length)  barriers,  it  involves  the  application  of  his  own  spherical 
wave  solution  (Thonasson,  1976)  to  the  ground  reflections  and  an 
Integration  of  the  total  field  over  the  face  of  the  barrier  itself7 
the  ground  reflections.  However,  in  order  to  obtain  a  solution, 
Thonasson  was  forced  to  Invoke  the  so-called  "Rayleigh"  or  "physical 
optics"  assumptions,  which  differ  only  slightly  from  the  Kirchoff 
assumptions. 

7  A  form  of  Bablnet's  Principle  is  needed  for  this;  see  Bouwkamp 
(1954,  pp.  51-52). 


RECEIVER  IMAGE 


Figure  1.4  Representation  of  the  Overall  Problem  as  a 

Superposition  of  Four  Diffraction  Paths  Using 
Source  and  Receiver  Images. 
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A  recent  contribution  by  Isel  (1980;  tee  also  Isei,  Embleton,  & 
Piercy,  1980)  is  very  similar  to  Jonasson's  treatment  except  that  it 

O 

allows  for  an  impedance  coating  on  the  faces  of  the  barrier. 

However,  the  Impedance  condition  is  handled  heurlstlcally  by  applying 
a  plane-wave  reflection  coefficient  to  the  second  term  in  the 
MacDonald  diffraction  solution  (the  term  ignored  by  Jonasson).  Hayek, 
et  al.  (1978),  on  the  other  hand,  used  Kendlg's  exact  diffraction 
solution  for  the  impedance-covered  half  plane,  merged  this  with 
Keller's  GTD,  but  applied  simple  plane-wave  theory  to  the  ground 
reflections.  The  present  study  has  extended  the  work  of  Hayek,  et  al. 
by  incorporating  into  their  model  an  original  solution  to  the  ground 
propagation  problem. 


Scope  and  Importance  of  Stud] 


This  presentation  will  begin  with  a  theoretical  treatment  of  the 
propagation  problem—that  is,  a  description  of  the  field  at  a  point 
(receiver)  in  the  half  space  above  an  infinite,  locally-reacting 
Impedance  plane  due  to  a  point  source  of  energy  located  arbitrarily  in 
that  half  space.  The  method  of  analysis  is  original,  although  an 
Important  step  early  in  the  derivation  is  similar  to,  and  was  inspired 
by,  a  step  Initially  used  by  Van  der  Pol  (1935)  in  his  electromagnetic 
work  and  later  adapted  by  Lawhead  and  Rudnick  (1951b)  to  the 
acoustical  case.  This  step,  which  involves  the  replacement  of  an 


fl 

Isei  apparently  forgot  to  account  for  one  image  term  in  his 
solution,  and  this  has  been  pointed  out  by  Fujiwara  (1980). 
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integrand  pole  by  an  equivalent  integral,  leads  to  a  further  integral 
which  ia  intractable.  Lawhead  and  Rudnick,  following  the  approach  of 
Norton  (1936,  1937),  made  a  major  approximation  so  that  the  integral 
could  be  solved.  The  approach  used  in  the  present  study  avoids  such 
an  approximation.  The  procedure,  which  will  be  described  in  detail  in 
chapter  II,  can  be  broadly  outlined  as  follows. 

First,  since  the  major  obstacle  with  the  integral  is  its 
complicated  exponential,  a  transformation  of  variables  is  made  which 
changes  the  exponential  into  a  more  useful  form.  The  denominator  of 
the  transformed  Integrand  is  then  expanded  in  a  Taylor  series  about 
the  origin  so  that  a  term-by-term  integration  can  be  performed.  Each 
Integral  in  the  resulting  series  is  then  solved  in  terms  of  parabolic 
cylinder  functions,  which  are  further  expressed  in  terms  of  the 
complementary  error  function.  The  final  formal  result  is  thus  an 

9 

asymptotic  series  (since  the  assumption  of  "large  kR"  was  necessary 
for  the  term-by-term  Integration)  containing  the  error  function  plus 
various  constants.  The  series  is  given  in  a  general  form,  including 
recursion  relations  for  calculating  the  coefficients.  Thus,  all  of 
the  higher-order  terms  in  the  asymptotic  series  can  be  written  out 
explicitly  and  their  values  can  be  calculated  readily. 

The  next  phase  in  the  study  couples  the  new  propagation  solution 
with  existing  half-plane  diffraction  theory  to  devise  a  practical 
model  for  handling  the  combined  barrier-ground  plane  problem.  The 

9 

Here,  k  is  the  familiar  propagation  constant  and  R  is  the 
distance  from  the  source  image  to  the  receiver. 
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development  of  this  model  Is  presented  in  Chapter  III. 

Nuaerlcsl  data  for  the  ground  propagation  problem  are  generated 
and  plotted  in  Chapter  IV,  Section  4.1.  Here,  the  sensitivity  of  the 
results  to  parameters  such  as  HkB,"  ground  impedance,  and  angle  of 
incidence  is  explored.  For  the  combined  barrier-ground  plane  problem, 
numerical  results  are  presented  in  Section  4.2,  where,  in  addition  to 
the  above  parameters,  the  dependence  on  barrier  height  and  diffraction 
angle  is  investigated. 

Although  the  major  scientific  contribution  of  this  research 
effort  is  the  asymptotic  series  solution  for  the  ground  propagation 
problem,  the  model  proposed  in  Chapter  III  for  the  combined  problem  is 
felt  to  be  an  important  engineering  contribution;  it  is  more  accurate 
than  either  Jonas sou's  or  Isei's  models  because  a  more  exact  ground 
propagation  theory  is  merged  with  a  more  exact  diffraction  theory. 
While  It  is  true  that  the  applicability  of  the  proposed  model  to  short 
barriers  is  questionable  (as  noted  by  Thomasson  in  regard  to 
Jonasson's  model),  a  numerical  integration  approach  such  as 
Thomasson's  (1978)— which  indeed  might  be  more  accurate  for  small 
barrier  heights— is  much  more  time  consuming  computationally,  and 
becomes  impractical  when  many  source-receiver  points  must  be 

i  ^  10 

analysed. 

^  Most  highway  noise  computer  prediction  schemes— a  common 
application  of  such  barrier-ground  plane  models— treat  the  flow  of 
vehicles  as  a  line  of  many  incoherent  point  sources.  As  several 
receivers  are  likely  to  be  included  in  any  one  scenario,  efficient 
computations  are  essential. 


CHAPTER  II 


MATHEMATICAL  FORMULATION  AND  SOLUTION 
TO  THE  PROBLEM  OF 

POINT  SOURCE  PROPAGATION  OVER  AN  IMPEDANCE  PLANE 

2.1  Overview 

The  exact  solution  to  the  problem  of  plane  wave  reflection  by  an 
impedance  plane  Is  rudimentary  and  well-known. 11  Until  recently, 
however,  it  was  customary  to  simply  take  the  results  of  the  plane  wave 
derivation  and  apply  them  to  the  point  source  (spherical  wave) 
problem.  Despite  the  lack  of  rigor,  the  solution  so  obtained  was 
consistently  confirmed  by  experiment,  except  when  the  sound  source  was 
very  close  to  the  reflecting  plane.  The  heuristic  derivation  of  this 
'‘plane-wave  solution"  to  the  point  source  problem  can  be  outlined  as 
follows . 

It  is  well-known  that  when  a  plane  wave  strikes  a  flat,  perfectly 

rigid  surface  at  an  angle  of  inclination  ^  ,  as  shown  in  Figure  2.1 

12 

(a),  the  resulting  field  at  a  point  (x,y)  can  be  written  as 

<j>(x,y)  .  ei[kxcostj>  -  kysinif/]  +  ei[kxcosi|<  +  kysin<|>]  (2a) 


11  See,  for  example,  Morse  and  Ingard  e  ^68,  p.  259). 

^  The  term  "field"  formally  refers  to  "velocity  potential,"  but 
since  harmonic  waves  are  most  likely  to  be  considered  here,  it  can 
be  taken  to  mean  "pressure"  as  well. 
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Here  k  is  the  wave  number,  or  propagation  constant,  and  k  ■  2uf/c, 
where  £  is  the  frequency  and  c  is  the  speed  of  sound  in  the  medium.  A 
common  factor  e  ,  where  io«2TTf  and  t  is  time,  has  been  ommitted  from 
this  and  the  following  equations.  The  first  term  in  Equation  (2.1)  is 
the  direct  contribution  to  the  field  from  the  source  acting  alone,  and 
the  second  term  represents  the  enhancement  of  the  field  due  to  the 
reflecting  plane.  The  total  field  at  the  point  (x,y)  would  be 
identically  expressed  if  the  plane  were  simply  replaced  by  an  "image" 
source,  as  indicated  in  the  figure. 

It  is  also  well-known  that  when  sound  energy  emanating  from  an 
ideal  point  source  S  strikes  a  rigid  plane  at  an  angle  ip,  as  shown  in 
Figure  2.1  (b),  the  resulting  field  at  the  point  P  is  given  by 


ikRL  ikR2 


where,  again,  the  first  term  is  the  direct  field  from  the  source  and 
the  second  term  can  be  thought  of  as  the  field  from  an  image  source. 

Now,  if  the  reflecting  plane  were  not  rigid  but  characterized  by 
a  complex  acoustic  impedance  Z,  the  plane  wave  solution  becomes 


$(x,y) 


gi[kxcosij;  -  kysimj/]  +  R  gi[kxcosij/  +  kysini^] 
S  P 


(2.3) 


Here  R  is  the 
P 


"plane  wave  reflection  coefficient"  given  by 


_  sinij*  -  8 
p  sini^  +  6 


(2.4) 
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in  which  3  “  pc/Z  is  the  normalized  acoustic  admittance  of  the 
13 

surface.  This  equation  is  exact  for  plane  waves;  it  results  from 
solving  the  full  boundary  value  problem.  The  solution  still  retains 
the  form  of  a  direct  term  plus  an  image  term,  but  now  the  "strength" 
of  the  image  is  no  longer  unity  but  R^. 

Since  Equation  (2.1)  transforms  into  Equation  (2.3),  then 
Equation  (2.2)  should  become 

IkR^  ikR2 

$(p)  -  ^-r—  +  R  (2.5) 

R1  P  R2 

for  the  spherical  wave  case  when  the  rigid  plane  is  replaced  by  an 
impedance  plane.  This  reasoning  has  been  employed  by  several 
researchers  to  support  the  use  of  Equation  (2.5)  to  solve  practical 
problems  involving  point  sources.  There  is  little  mathematical  basis 
for  this  equation,  and  it  is  not  surprising,  therefore,  that  when 
acoustic  energy  from  a  real  source  ceases  to  resemble  "plane" 
waves — at  grazing  incidence,  for  example— the  solution  falls. 

What  is  needed  is  a  "spherical  wave  reflection  coefficient”  for 
the  point  source  problem.  That  is,  the  solution  should  take  the  form 


ikRx  ikR2 
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The  normalization  is  by  the  characteristic  Impedance  Pc  of  the 
medium,  where  p  is  the  density.  The  normalized  acoustic  impedance, 
sometimes  called  the  "specific  acoustic  impedance,"  is  Z„  -  Z/Pc  and 
hence  6  -  1/Z„  .  tJ 

N 


where  Q  is  derived  analytically  from  the  formal  boundary  value 
problem.  Alternatively,  a  solution  may  be  sought  in  the  form 
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$<P) 


ikR  ikR2 


Correction 

term 


(2.7) 


where  the  correction  term  approaches  zero  in  the  limit  when  the  source 
is  very  far  from  the  boundary  but  becomes  significant  when  the  source 
is  close  to  the  impedance  plane.  Both  of  the  above  forms  of  solution 
will  be  derived  in  this  study,  but  the  practical  selection  of  one  form 
over  the  other  will  depend  on  the  intended  application. 

Perhaps  some  physical  insight  about  the  problem  can  be  gained 
before  the  theoretical  treatment  begins.  Figure  2.2  (a)  and  (b) 
depicts  three  elemental  segments  along  a  single  wavefront  of  either  a 
plane  wave  or  a  spherical  wave,  respectively.  The  rays  perpendicular 
to  the  two  wavefronts  follow  the  progression  of  the  segments  A,  B,  and 
C  to  their  "striking  points”  along  the  boundary.  A',  B',  and  C'.  Of 
course,  each  segment  on  the  wavefront  strikes  the  boundary  at 
different  times,  but  this  is  not  important  here.  Rather,  the  angle  of 
incidence  with  which  they  strike  the  boundary  is  the  important  factor. 
For  the  plane  wave,  each  elemental  portion  of  the  wavefront  hits  the 
boundary  at  the  same  angle,  but  for  the  spherical  wave,  this  is 
clearly  not  the  case.  With  the  assumption  of  local  reactance,  each 
ray  "sees”  a  surface  impedance  proportional  to  its  normal  component  of 
particle  velocity  only.  This  component,  therefore,  has  the  same  value 
for  each  ray  in  Figure  2.2  (a)  but  different  values  for  each  ray  in 


Segments  of  a  Wavef 
erical  Wave. 
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Figure  2.2  (b).  Consequently,  the  phase  change  upon  reflection  is 
the  same  at  the  points  A",  B' ,  and  C'  for  the  plane  wave  but  varies 
from  A'  to  E'  to  C"  for  the  spherical  wave.  The  plane  wave  is 
reflected  as  a  plane  wave,  but  the  spherical  wave  loses  its  symmetry. 
When  the  ground  plane  is  perfectly  rigid  (or  pressure  release),  the 
angle  of  incidence  is  immaterial  since  the  phase  change  is  always  180 
degrees.  Hence,  for  a  rigid  plane  (or  pressure  release  plane)  only, 
the  spherical  symmetry  is  preserved,  and  the  image  method  is  ex*  .t. 

2.2  Formulation  of  the  Problem 

An  Idealized  point  source  is  located  at  a  fixed  but  arbitrary 
position  in  the  half  space  above  a  perfectly  flat  ground  plane  that 
extends  to  infinity  in  its  two  dimensions.  The  source  radiates 
acoustic  energy  into  the  uniform,  quiescent  medium  that  fills  the 
space  above  the  plane.  The  resulting  field  at  a  point  receiver,  also 
located  at  an  arbitrary  position  above  the  plane,  is  desired.  By 
definition,  the  ground  plane  is  "locally  reacting;"  that  is,  the  value 
of  the  acoustic  impedance  of  the  surface  facing  the  source  is 
independent  of  the  angle  of  Incident  energy  upon  the  plane,  and  it 
may,  therefore,  be  fully  characterized  by  its  value  at  normal 
incidence.  Also,  in  this  ideal  case,  the  proximity  of  the  ground  has 
no  effect  on  the  radiation  or  "strength"  of  the  source. 

Mathematically,  the  problem  stated  above  can  be  described  by  an 
inhomogeneous  partial  differential  ("wave")  equation  subject  to 
homogeneous,  mixed  boundary  conditions.  Allowing  only  single 
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frequency,  steady-state  excitations  will  reduce  the  governing 
differential  equation  to  the  Helmholtz  equation.  The  general  solution 
to  such  a  boundary  value  problem  comprises  a  particular  solution  to 
the  inhomogeneous  differential  equation  plus  a  separate  solution, 
containing  unknown  coefficients,  to  the  homogeneous  differential 
equation.  The  total  solution  must  satisfy  the  boundary  conditions. 

For  this  reason,  it  is  advantageous  to  express  the  particular  and 
homogeneous  parts  of  the  solution  in  coordinates  that  are  easily 
matched  to  the  geometry  of  the  bounding  surface  itself. 

2.2.1  The  Particular  Solution:  A  Point  Source  in  a  Free  Field 

In  the  absence  of  the  bounding  impedance  plane,  the  problem 
reduces  to  a  point  source  of  sound  radiating  into  an  unbounded  medium. 
This  problem  is  most  easily  described  in  terms  of  spherical 
coordinates,  and  the  solution  is  derived  in  many  standard  textbooks. 

It  is  simply: 

lkR 

*<P)  -  ,  (2.8) 

where  R  is  the  distance  between  the  point  source  and  an  arbitrary 
field  point  P.  Although  this  uncomplicated  form  of  solution  is 
appealing,  it  can  not  be  matched  to  the  planar  boundary  conditions. 

2. 2.1.1  The  Field  from  a  Point  Source  in  Cartesian  Coordinates. 
The  boundary  plane  can  be  specified  in  Cartesian  (x,y,z)  coordinates 
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(for  example,  z-0) ,  and  therefore  a  solution  to  the  Helmholtz  equation 
in  this  system  would  be  more  useful.  The  Inhomogeneous  equation  is 


V2*(x,y,z)  +  k24>(x,y,z)  -  -4tt  6(x)6(y)6(z-zQ)  ,  (2.9) 


is  defined  as  the  Fourier  transform  of  4»(x,y,z),  and 

QO 

*(x,y,z) - ~  fff  #<e,Tuc>  ei(Cx  +  ny  +  ?z)d?dndC  (2.12) 

(2ir)3  >11 

is  then  defined  as  the  Inverse  Fourier  transform. 

Now,  applying  the  transform  to  Equation  (2.9)  and  employing  the 
following  properties  of  the  Fourier  integral  and  the  delta  function, 
respectively: 


Figure  2.3  The  Cartesian  Space  Geometry,  (a)  Source  and 
Receiver  Coordinates,  (b)  Propagation  Vector 
for  a  Single  Elemental  Plane  Wave  with  Arbitrary 
Direction  in  the  Space. 
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dt  -  -iyf(Y) 


(2.13) 


|  6(t-tQ)  f (t)dt  -  f(tQ)  ,  (2.14) 

the  field  in  transform  space  becomes 

"i?Zo 

$(C,n,0  -  -—--f - - - —  (2.15) 

(r  +  n  +  r  -  o 

Applying  the  inverse  transformation,  Equation  (2.12),  to  this  will 
yield  the  acoustic  field  in  the  original  coordinate  system  due  to  the 
point  source  at  (0,0,z0).  Thus, 

i  rTr  e1^* +  w  +  ?(z“2o^ 

*(x,y,z)  - r  - r - = - = - 5 — d£dnd£  .  (2.16) 

(2")Z  Jjj  (C  +  n  +  C  -  kZ) 

The  parameters  £,  n»  and  C  may  be  thought  of  as  components  of  a 
propagation  vector,  in  which  case  the  integrand  in  Equation  (2.16) 
takes  on  the  form  of  a  plane  wave.  More  precisely,  since  £,  n,  and  £ 
may  have  arbitrary  and  independent  values,  the  integral  describes  a 
superposition  of  an  infinite  number  of  plane  waves  of  different 
amplitudes  propagating  in  different  directions.  Equation  (2.16)  can 
be  written  in  a  form  that  portrays  these  ideas  more  clearly: 

1  fff 

$(R)  -  $(x,y,z)  -  - j  III  ~2 - 2  dWnd?  (2.17) 

(2ir)  _0o  X  —  k 
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rs 

where  Che  propagation  vector  X  -  £ig  +  ni^  +  has  a  magnitude 

2222  * 

X  ■  £  +  n  +  5  ,  and  the  radius  vector  R  ■  xi  +yi  +  (z-z  )i  has  a 

x  y  02 

2  2  2  2 

magnitude  R  *  x  +  Y  +  ( Z~ZQ )  •  A  single  component  plane  wave 

with  propagation  vector  X  is  shown  in  Figure  2.3  (b)  to  illustrate 
these  concepts  geometrically.  Finally,  it  should  be  noted  that,  with 
a  few  simple  transformations,  Equation  (2.16)  can  be  integrated  (for 
one  treatment,  see  Stratton,  1951,  pp.  577-582)  to  generate  the 
familiar  point  source  representation  of  Equation  (2.8). 


2*2. 1.2  The  Field  from  a  Point  Source  in  Cylindrical 
Coordinates .  The  boundary  plane  could  be  specified  by  a  single 
coordinate  surface  in  cylindrical  (r,0,z)  as  well  as  Cartesian 
coordinates  (for  example,  z«0),  and  therefore  a  solution  to  the 
Helmholtz  equation  in  cylindrical  coordinates  could  be  conveniently 
applied  to  the  boundary  conditions.  The  Helmholtz  equation  for  a 
source  located  at  a  point  z*zQ  along  the  cylindrical  axis  takes  the 


V  4(r,z)  +  k  $(r,z)  - 


-2  6(r)  6 (z-z  ) 


(2.18) 


where  the  Laplacian  operator  in  cylindrical  coordinates  is 

c2  11  ,  3  x  .  32 

V  "  7  3?  (r  37J  +  rr 

+  (2a9) 
3r2  +  r3r  +  a,2 

Because  the  source  has  been  located  directly  on  the  z-axls  of  the 
coordinate  system,  there  is  no  azimuthal  angle  dependence  in  either  of 
the  preceding  equations. 


d 


Nov,  applying  the  one-dimensional  Fourier  transfora  to  the  axial 
distance  z,  where 
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Aa.z) 


-v(z-zo) 


;  (z-zq)  >  0 


+v(z-z) 

o 


;  (z-zq)  <  0 


(2.30) 


2  i  !« 

which,  under  the  constraint  Re(V)  ■  Re[(A  -k  )  j  >0,  can  be  written 


♦*(A,z)  -  ■£■ 


-V  z-z. 


(2.31) 


Finally,  the  inverse  Hankel  transform.  Equation  (2.21),  is  applied  to 
Equation  (2.31)  to  obtain  the  desired  field  at  the  point  (r,z)  in 
cylindrical  coordinate  space: 


w  — 

lkR  r 


-V  Z-Z. 


- - XjQ(Xr)  dA  , 


(2.32) 


where  R  -  /r2  +  (z-zQ)2  . 

This  result  is  the  particular  solution  to  the  inhomogeneous 
Helmholtz  equation.  In  other  words.  Equation  (2.32)  represents  a 
point  source  radiating  into  a  free  field.  From  fundamental 
principles,  it  is  known  that  functions  of  the  form 

ein®  JQ[(k2-h2)V]  eillz  ,  where  h  and  k  are  propagation  constants  and  9 
is  the  azimuthal  angle,  represent  elementary  cylindrical  waves,  which 
are  solutions  to  the  homogeneous  Helmholtz  equation.  The  point  source 


The  wavenumber  k  is  assumed  to  be  complex  in  general;  that  is, 
the  medium  is  assumed  to  possess  a  small,  but  finite,  amount  of 
dissipation. 


I 
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In  Equation  (2.32)  thus  represents  an  infinite,  weighted  sum  of  these 
elementary  cylindrical  waves  propagating  symmetrically  about  the  z- 
axis  (n-0).  Sommerf eld  (1909)  was  the  first  to  derive  this  result, 
and  the  representation  of  Equation  (2.32)  is  sometimes  called  the 
"Sommerf eld  point  source.” 


2.2.2  Total  Solution:  Effects  of  the  Boundai 


Each  of  the  representations  of  a  point  source  in  free 

space— equation  (2.16)  in  terms  of  plane  waves  and  Equation  (2.32)  in 

terms  of  cylindrical  waves — is  a  solution  to  the  non— homogeneous 

Helmholtz  equation.  When  a  boundary  such  as  the  impedance  ground 

plane  is  present,  a  second  term  must  be  added  to  the  free-space 

solution  to  account  for  reflections.  The  total  field  can  be  written 

as  $tot  ■  +  $r  ,  where  the  components  may  be  termed  the  "incident 

field"  and  the  "reflected  field,"  respectively.  Because  the 

lkR 

particular  solution—  *  e  /R  ,  or  either  of  its  elementary  wave 
expansions— is  unique,  the  reflected  term  can  not  also  be  a 
solution  to  the  non-homogeneous  equation,  but  instead  satisfies  the 
homogeneous  Helmholtz  equation 


V2$(r,z)  +  k2$r(r,z) 


(2.33) 


As  such,  the  solution  for  the  reflected  field  will  contain 
arbitrary  constants  which  can  be  determined  by  applying  the  boundary 


conditions 
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2. 2. 2.1  Total  Integral  Solution  via  Elementary  Wave  Functions. 
The  free-space  solutions  developed  in  the  previous  sections  also 
satisfy  the  the  homogeneous  Equation  (2.33)  at  all  points  except  at 
the  source  itself  (R»0).  The  elementary  plane  and  cylindrical  waves 
are  therefore  eigenfunctions  of  the  problem,  and  the  reflected  field 
can  be  constructed  from  an  infinite  summation  (integral)  of  these 
eigenfunctions,  with  arbitrary  weighting  functions  to  be  determined 
from  the  boundary  conditions. 

The  expansion  of  the  reflected  field  in  terms  of  elementary  plane 
waves  was  the  approach  used  by  Weyl  (1919)  for  the  electromagnetic 
problem  and  later  adopted  by  Ingard  (1951)  for  the  acoustic  case.  By 
defining  y  ■  [£2  +  r\^  -  with  Re  Y  >  0,  Equation  (2.16)  can  be 

transformed  into 


4>(x,y,z) 


i(£x  +  ny) 


d£dn 


i£(z-zo) 

e _ 

(C  +  iY)(C-iY) 


dC  .  (2.34) 


Using  residue  methods  for  the  5-integration,  this  can  be  reduced  to 


$(x,y,z) 


00  -y  |  z-z  | 

_L  ff  e _ i(£x  +  ny) 

2tt  jj  y  e 

—CO 


d^dn 


(2.35) 


In  this  form,  the  elementary  plane  waves  are  decaying  exponentially 
with  z,  and  Weyl  termed  these  "inhomogeneous”  plane  waves.  Noting 
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that  C  +0  +  (iY)  ■  k  ,  the  propagation  constant  It  can  be  written 

as  a  vector  with  components  k  *  £*  +  n!  +  (iy)I  •  This  suggests 

x  y  z 

the  transformations  iy  ■  kcosa,  £  -  (ksina)cosB,  and  n  “  (ksina)sing. 
The  resulting  field  at  the  point  (x,y,z),  then,  from  the  point  source 
located  at  z-z^  can  be  written  as 

<$(x,y,z)  - 


Tr/2-i00  2n 


ik 

27T 


ik[xsinacosS + yslnasinB +  |z-z  |cosa] 


sinadadB 


(2.36) 


The  final  step  is  to  assume  that  the  reflected  field  "originates"  from 
the  image  point  at  z  «  -zQ  and  thus 


<frr(x,y,z)  - 

ir/2-i00  2tt 
ik  f  .  .  . 

15  j  j  V"!' 


ik[xsinacos0  +  ysinasin6+ jz+z  |cosa]  (2.37) 

L  •  sinadadB 


where  the  unknown  function  fr(ot)  must  be  determined  from  the  boundary 
conditions.  For  the  present  case  of  an  impedance  boundary,  f^(a)  is 
simply  the  plane  wave  reflection  coefficient  given  by  Equation  (2.4). 

The  expansion  of  the  reflected  field  in  terms  of  elementary 
cylindrical  waves  was  the  approach  used  by  Sommerfeld  (1909)  and  Van 
der  Pol  (1935)  for  the  electromagnetic  problem,  and  later  adapted  by 
Rudnick  (1947)  to  the  acoustical  case.  To  express  the  reflected  field 
in  integral  form,  the  results  of  Section  2. 2. 1.2  are  used, 
specifically  Equation  (2.32).  Assuming  the  reflected  field  originates 
from  an  image  point  source  at  z— zQ,  the  field  can  thus  be  written  as 
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oa  -V(Z-Z  ) 

<&(r,z)  -  |  gr(X)  - - - -  XJo(Ar)  dA .  (2.38) 

o 

where  the  unknown  function  g^(V)  oust  be  determined  from  the  boundary 
conditions. 

The  boundary  conditions  for  a  locally  reacting  surface  is  that 
the  ratio  of  the  pressure  p  to  the  normal  component  of  the  velocity  vq 
equals  the  acoustic  impedance  of  the  surface.  That  is, 

.  Z  .  (2.39) 

v  (r,z  v  (r,o) 

n  n 

z-0 

Since  the  pressure  is  related  to  the  velocity  potential  by  the 

fundamental  equation  p  -  p3$/3t,  or  p  -  -imp*  for  harmonic  waves,  and 

since  v  *  -3$/3z,  the  boundary  condition  can  be  expressed  in  the  form 
n 

3<t>.(r >.9I  +  lk04 (r , 0)  -  0  .  (2.40) 


Here  0*pc/Z  is  the  normalized  acoustic  admittance  of  the  surface. 

When  the  total  solution  ♦  -  ^  +  *r  is  substituted  into 

Equation  (2.40),  the  unknown  function  g^( A)  can  be  determined  to  be 


grtt) 


v  +  ikS 
v  -  ikS  ‘ 


(2.41) 


Therefore,  the  integral  form  of  the  total  solution  is 


$(r,z) 


00 


O 


V  +  ikg 

V  -  ikg_ 


-v(z+zQ) 

e _ 

V 


XJ  (Xr)  dX 
o 


(2.42) 


2. 2. 2. 2  Total  Integral  Solution  via  Hankel  Transforms.  The 
result  just  obtained  can  be  derived  In  a  more  stralghforvard  manner 
using  Hankel  transforms.  The  total  solution  to  the  transformed 
Helmholtz  equation  (2.24)  comprises  the  particular  solution. 

Equation  (2.30),  plus  a  homogeneous  solution  containing  two  arbitrary 
constants  to  be  determined  from  two  boundary  conditions.  Thus,  the 
field  in  transform  space  takes  the  form 

-v|z-zQ| 

$*(X,z)  -  A(X)  e-VZ  +  B(X)  eVZ  +  * - - -  .  (2.43) 

The  first  boundary  condition  that  must  be  met  is  the  Sommerfeld 
radiation  condition  specifying  that  the  field  must  vanish  as  z  . 
Because  the  constraint  Re(v)>0  has  already  been  imposed,  the  constant 
B  must  be  set  equal  to  zero.  The  second  boundary  condition  is  that 
dictated  by  the  local  reaction  assumption,  namely.  Equation  (2.40). 
Applying  the  Hankel  transform.  Equation  (2.20)  to  the  boundary 
condition  yields 

^-4^1+  ikB$*(X,0)  -  0  .  (2.44) 

aZ 

* 

Substituting  $(X,z)  from  Equation  (2.43)  into  Equation  (2.44)  and 
writing  the  absolute  value  quantity  as  |z-Zq|  *  (z^-z)  for  z-»-0  yields 


the  desired  value  of  the  constant: 
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(2.45) 


The  solution  for  the  field  in  transform  space,  then,  becomes 


♦  (X,z) 


-v(z+zq) 

e 

v 


+ 


V 


(2.46) 


Finally,  taking  the  inverse  Hankel  transform.  Equation  (2.21),  of  this 
expression,  yields  the  integral  form  of  the  total  solution.  Hence, 


oo  — v  I  z-z  I  •  — v (z+z  ) 

#  I  0  1  /  ■■  ■_  Q 

dX  +  j  *— - - XJo(Xr)  dX 

o  o 


oo 


O 


-v(r**zo) 

1 - - - XjQ(Xr)  dX 


(2.47) 


which  is  identical  to  the  previous  results  of  Equation  (2.42) 

2.3  Derivation  of  the  Asymptotic  Series  Solution 

The  integral  in  the  exact  solution,  Equation  (2.47),  appears  to 
be  intractable  due  to  the  presence  of  a  simple  pole  at  v  - 
(X^-k^)*5  ■  +ikg  and  branch  points  at  X-  ±k.  Virtually  all  of  the 
researchers  who  have  addressed  the  spherical  wave  propagation  problem 
over  the  past  century  have  been  faced  with  this,  or  a  very  similar, 
integral,  and  all  have  necessarily  resorted  to  approximation  methods 
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for  its  evaluation.  Those  methods  based  on  an  asymptotic  approach,  in 
which  a  critical  parameter  such  as  "kR"  is  assumed  to  be  large,  have 
proven  to  be  the  most  accurate  and  to  provide  the  most  useful  results. 
In  these,  the  real  axis  path  of  Integration  is  continued  into  a 
contour  in  the  complex  plane  that  invariably  must  be  deformed  around 
the  branch  point  and  branch  cut.  It  is  to  the  resuting  branch  line 
Integral  that  the  asymptotic  technique—usually  the  "method  of 
steepest  descents" — is  applied.  In  addition,  proper  attention  is  paid 
to  the  simple  pole,  which  may  or  may  not  be  enclosed  by  the  final 
contour  depending  on  the  specific  nature  of  the  problem. 

The  approach  in  this  study  is  different.  The  pole  is  replaced  at 
the  outset  by  an  integral  representation,  and  the  branch  line 
integration  is  avoided  through  an  appropriate  transformation  of 
variables.  However,  the  "asymptotic"  assumption  is  still  necessary, 
but  it  is  invoked  later  in  the  development  when  the  integral  is 
expanded  into  a  Taylor  series.  The  final  form  of  solution  will  be  an 
asymptotic  series,  for  which  the  coefficients  may  be  calculated  by 
recursion  formulae. 

2.3.1  Replacement  of  the  Pole 

By  re-writing  (v+y)/(v-y)  ■  1  +  2Y/(v-y),  the  final  solution 
given  in  Equation  (2.47)  for  the  field  at  the  point  (r,z)  can  be 
expressed  as 
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ikR1  . ikR2 


00  -v(z+z  ) 

(  .  o 


t(t,z)  -  +  S1—  *  2y  |  '(v_Y)v  Wo(Ar)  dA 


(2.48) 


in  which  Y"ik8,  R,  is  the  distance  from  the  source  at  (0,Zq)  to  the 
point  (r,z),  and  R^  is  the  distance  from  the  image  point  at  (0,-Zq)  to 
the  point  (r,z).  The  geometry  of  the  problem  is  illustrated  in  Figure 


From  the  fact  that  1/p 


f 


dt  for  Re(p)>0,  the  pole  in 


the  Integrand  of  Equation  (2.48)  can  be  transformed  asJ 


w 

■hr  "I 


e  (v  -  y)R|  .  Re  [V  -  Y]  >  0 


(2.49) 


However,  it  turns  out  that  the  constraint  Re(v-Y)>0  is  not  very 
16 

convenient,  and  a  better  transformation  would  be 


(v  -  r) 


J  e"i(V  "  Y)n  dn  ;  Re  [i(v  -  Y)1  >0  (2.50) 


so  that 


w 

— y-  i  |  e"i(v  "  Y)n  dn  ;  Re  [i(v  -  y)J  0  .  (2.51) 


13  This  step  is  based  on  a  similar  one  first  used  by  Van  der  Pol 
(1935)  and  later  adopted  by  Lawhead  and  Rudnick  (1951). 

**  The  constraint  translates  into  Re(v)  +  klm(g)  >  0.  Although 
Re(v)>0,  the  imaginary  part  of  g  may  be  positive  or  negative,  and  the 
condition  Re(v~Y)>0  can  not  be  guaranteed. 
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Figure  2.4  Geometry  of  the  Propagation  Problem  in 

Cylindrical  Coordinates.  Source-Receiver 
Distance  R.,  Image-Receiver  Distance  Rg, 
Reflection  angle  \ \i,  and  Ground  Admittance  8. 
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Writing  v  -  vr  +  iv£  and  Y  -  ik(Bf  +  iB^  -  -kB1  +  ikBr,  where  the 
subscripts  "r"  and  "i"  indicate  "real"  or  "imaginary"  parts,' 
respectively,  then  the  constraint  Re[i(v-y)]  >  0  translates  into 

(“Y^  +  kBr)  >  0.  Because  the  real  part  of  the  impedance  Z«R+iX  for 

2  2 

physical  materials  must  be  positive,  then  6^  -  R/(R  +X  )  must  also  be 

positive.  Thus,  in  order  to  satisfy  the  constraint  in  Equation 

(2.51),  it  is  sufficient  that  be  negative.  That  this  is  Indeed  the 

case  requires  a  closer  look  at  the  parameter  v. 

Assuming  that  the  medium  is  not  Ideally  lossless,  the  wavenumber 

k  may  be  written  as  a  complex  quantity  k  ■  kr  +  ik^  and  because 
ikR 

e  /R  must  decay  for  large  R,  the  imaginary  part  k.  must  be  positive. 
Furthermore,  if  the  symbol  "/“"  is  adapted  here  to  indicate  the 
positive  real  square  root  (often  called  the  "principal  root"),  then 
v  •  A.  -  k2  -  /"(^-k^+k2)  -  i(2kfki) .  The  real  part  of  the 
quantity  under  the  radical  can  be  positive  or  negative,  but  the 
imaginary  part  will  always  be  negative.  Hence,  the  root  must  lie  in 
the  fourth  quadrant;  that  is,  v^O.1^ 

The  replacement  of  the  pole  via  Equation  (2.51)  is  thus  valid  for 
general  Impedances  (springlike  or  masslike),  and  when  the  result  is 
substituted  into  Equation  (2.48),  the  integral  term  becomes 


This  can  also  be  shown  by  considering  the  assumed  time 
dependence .  Thus , 

#v|*-20|e-iu>t  .  #-vr|z-*0le(-vi|z-z0|-u)t)  . 

therefore,  to  ensure  outward  propagation,  v^O. 


w  w 

!1  ■  “V  {{ 


-v(hfin)  I  lYn 

- - - -  XjQ(Xr)  dX(  e  Tn  dn 


i 


(2.52) 


o  o 


where  h^Cz+Zg).  The  term  In  brackets  resembles  the  Sommerfeld 
representation  for  a  “point  source,"  Equation  (2.32),  and  therefore 


[i  - 2i^  ]  [t^-]  'lTn 


dn 


(2.53) 


where 


R' 


/?  +  (h+in)2  -  A*  +  h2  +  2ihn  -  n2 


(2. 54) 


2  2  2 

From  Figure  2.4  it  can  be  seen  that  h/R^  »  sin^  and  r  +h  -R2,  thus 


R'  -  R„ 


1  +  2*n  sinij;  rj_ 
“2  R2 


(2.55) 


where  Re(R“’)  >  0. 

The  integral  in  Equation  (2.53)  now  takes  the  form 


ikR 

*  e 

—  f 

1  +  2in-?inip  _  n 
.1  R2  R2 

+  iyn 

J 

n 

.  2in  sin $  & 

* 

%2  n2 

z  R2J 


dn 


(2.56) 


or 
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t1-  -21c0  j 


»  ikR2f(l  -  52  +  2i£  sin^)*5  +  IBS] 


2  ST 

(1-5  +  2iC  sintj;)^ 


d5 


(2.57) 


where  the  transformation  5-n/R2  has  been  applied.  Although  the  pole 
has  been  "removed”  from  1^,  the  Integral  is  still  cumbersome,  and  a 
transformation  of  variables  will  now  be  made  in  order  to  simplify  the 
exponential  term. 


2.3.2  Variable  Transformation 

The  functional  form  of  a  single  variable  transformation  5 “€(t)  on 
a  contour  integral  is 


f<5>  dC 


J  f[£(t)]^^-dt  , 


(2.58) 


or,  more  specific  to  the  case  at  hand, 

|  «a£(5)  F(5)  dC  *  j  •*ft£le>]  at  .  (2.59) 

ce  Ct 

where  and  are  arbitrary  contours  in  complex  5-space  or  complex 
t-space,  respectively. 

The  exponential  in  the  field  solution.  Equation  (2.57),  can  be 


put  into  a  simpler  form  by  the  variable  transformation 
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A  -  52  +  2i£  sinip  +  IBS  -  -(t2  +  2Bt  +  C)  .  (2.60) 

With  this  transformation,  the  original  real  axis  path  0<£<«will  in 
general  become  a  contour  in  the  complex  t-plane. 

To  determine  the  values  of  the  coefficients  B  and  C,  two 
constraints  must  be  placed  on  the  transformation.  The  first  is  that 
the  transformed  contour  should  start  at  the  origin;  that  is,  t(5*0)*0. 
This  condition  Immediately  fixes  the  value  O-i  (again,  the  root  on 
the  left  hand  side  of  Equation  [2.60]  is  taken  so  that  the  real  part 
is  positive).  The  second  constraint  on  the  transformation  is  less 
trivial.  Of  course,  if  the  transformation  were  matched  at  any  other 
point  t(£g)-tg,  the  unknown  coefficient  B  could  be  determined. 

However,  the  most  judicious  choice  for  this  matching  point  is  the  so- 
called  "saddle  point"  of  the  original  Integrand. 

2. 3. 2.1  Matching  Transformation  at  the  Saddle  Point.  The  point 
at  which  the  derivative  f '(5)*0  is  called  the  "saddle  point"  and  is  a 
fundamentally  important  quantity  in  the  asymptotic  technique  of 
"saddle  point  Integration,"  also  called  the  "method  of  steepest 
descent."  At  the  saddle  point,  the  integrand  takes  its  maximum  value 
and  falls  off  rapidly  to  either  side  in  the  directions  of  "steepest 
descent."  The  power  of  this  method  lies  in  the  fact  that  if  the 
integration  is  performed  along  the  path  of  steepest  descent,  most  of 
the  contribution  to  the  integral  comes  from  the  immediate  vicinity  of 
the  saddle  point. 

The  saddle  point  of  the  integrand  in  Equation  (2.57)  satisfies 


51 


[i(l  -  52  +  2i§  sliwj#)5*  +  1651  -  0  ,  (2.61) 

and  the  saddle  point  tg  for  Equation  (2.57)  satisfies 

-f-  i[-t2  -  2Bt  +  1]  -  0  ,  (2.62) 

dt 

where  the  parameter  a  In  Equation  (2.59)  has  been  set  equal  to  kRj* 
Therefore,  solving  the  previous  two  equations,  the  saddle  points  are 
found  to  be 


(2.65) 


Note  that  at  this  point  in  the  derivation  there  is  nothing  to  govern 
the  selection  of  the  root  in  Equation  (2.65). 


2. 3. 2. 2  The  Transformation  of  the  Integral.  The  transformation 


of  the  exponential  factor  is  complete;  that  is. 


ikR.  [  (1  -  £2  +  2i£  sintii)5*  +  i0£]  -ikR,[t2  +  2Bt  -  1] 


(2.66) 


With  reference  to  Equation  (2.59),  the  next  step  in  the  transformation 
involves  the  function  F(£),  which  is  defined  here  as 


(2.67) 


4  -  5  +  2iC  sintp 


In  order  to  find  F[£(t)]  where 


F[C(t)]  - 


(2.68) 


4  -  C  (t)  +  2i£(t)  slnip 


the  quantity  £  must  be  expressed  as  a  function  of  t.  In  a 
straightforward  manner,  Equation  (2.60)  can  be  solved  for  £(t)  as 


£(t)  .  - i-—  {$(-t2  -  2Bt  +  1)  +  sin  if 

(1  -  0  ) 


(2.69) 


±[(-t2  -  2Bt  +  1  +  sin^)2  -  /l  -  02^  cost//]15} 


The  quantity  in  the  square  bracket  can  be  factored  so  that 
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where 


5(tX - ^-5-  {ec-t2  +  2Bt  +  1)  +  sini^ 

(1  -  8  )  1 


f8  +  simjj]  , 

t2  2Bt\ 

L  B  J 

H  H  J 

(2.70) 


H 


1  +  8  sint|>  + 


(2.71) 


and  where  the  sign  of  the  root  in  Equation  (2.69)  had  been  selected  to 
make  £(t*0)»0,  as  required.  As  a  check  on  the  transformation  in 
Equation  (2.70),  it  is  easy  to  confirm  that  C(tQ)*?Q.  Now  when  £(t) 
from  Equation  (2.70)  is  substituted  into  Equation  (2.68),  the 
integrand  function  F[£(t)]  can  be  determined.  Hence, 


F[C(t)] 


(1  -  (ft _ 

[a(t)  +  8  sin*]  -  B(t  +  B)  /yTt7 


(2.72) 


in  which  the  following  shorthand  notations  have  been  introduced: 


a(t) 


2Bt  +  1 


Y(t) 


(2.73) 


The  final  step  in  the  transformation  of  the  integral  is  to  derive 
an  expression  for  d£(t)/dt.  This  is  effected  by  differentiating  the 
expression  for  ?(t)  in  Equation  (2.70),  and  therefore 


d£(t)  _  i 

• 

-8(t  +  B) 

8  +  siniiJ 

B 

/y(t)'  +  [a(t)  +  8  sin^] 

dt  (1-32) 

k 

6  +  sini)) 
B 

/y(t)' 

(2.74) 


The  product  F[£(t)]d£(t)/dt,  then,  results  in  the  relatively  simple 
expression 


FlCCt)]  - £  IB  -  .  (2.75) 

(3  +  sin<|>)  /l  -  ^  ^ 

With  the  transformation  thus  complete,  the  original  integral  of 
Equation  (2.57)  can  be  written 


I 


1 


ikR 

-4ik8Be 


(8  +  simjj) 


-ikR  [t2  +  2Bt] 
(  & 


e 


dt 


(2.76) 


with  H  defined  in  Equation  (2.71)  and  B  in  Equation  (2.65).  The 
integral  is  taken  over  the  path  Ct,  the  transformed  contour  in  the  t- 
plane.  The  root  in  the  expression  for  B  will  be  fixed  by  the  choice 
of  this  integration  path. 

2. 3. 2. 3  The  Path  of  Integration  in  the  t-plane.  As  the  value  of 


£  goes  from  0  to  +00 along  the  real  axis,  the  transformed  variable  t 
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will  trace  out  a  path  in  the  complex  t-plane.  The  transformation. 
Equation  (2.60),  can  be  readily  solved  for  t(£),  which  will  describe 
this  path;  thus. 


t<€)  -  -B  ±  [B2  -  Vtt)]* 


(2.77) 


where 


V(5)  -  -1  +  A  -  £2  +  2isin +  i0£  .  (2.78) 

The  constraint  t(£-0)  -  0,  or 

t(0)  ■  -B  ±  [B2]*5  (2.79) 


dictates  the  selection  of  the  root  in  Equation  (2.77). 

The  parameter  B  =  ±i(l  +  0sini{/  -  /i-^'costy)*5  can  be  neglected  in  ' 
comparison  with  V(£)  for  where  R  is  some  large,  positive 

constant.  Thus,  from  Equation  (2.77), 


t(5  -  R)J  -  ±[1  -  A  -  R2  +  2iRsinip  -  10R]1* 
R-x» 

-  ±[-  -  i(0r  +  i0i)  R]*5 

-  ±[8^  -  i(R  +  erR) j1* 


(2.80) 


The  imaginary  part  of  the  quantity  inside  the  brackets  is  always 
negative,  but  the  real  part  can  be  negative  or  positive.  Thus,  if 
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$j<0,  the  quantity  in  brackets  is  confined  to  the  third  quadrant,  and 
the  two  roots  will  lie  in  the  seventh  and  third  octants.  If  B^>0,  the 
bracketed  term  is  in  the  fourth  quadrant,  and  therefore  the  two  roots 
will  lie  in  the  eighth  and  fourth  octants.  These  relationships  are 
illustrated  in  Figure  2.5  (a). 

The  octants  that  have  been  identified  represent  potential  closure 
regions  for  the  integration  contour  in  the  complex  t-plane.  That  is, 
the  original  path  along  the  real  axis  will  connect  via  a  circular  arc 
of  large  radius  to  the  point  t(R)  expressed  in  Equation  (2.80). 
Refering  to  Figure  2.5  (a),  it  is  apparent  that  the  most  convenient 
choice  for  closure  is  the  eigth  octant  for  6^>0  and  the  seventh  octant 
for  Bj^O.  In  either  case,  the  final  closed  contour  will  enclose  a 
region  entirely  within  the  fourth  quadrant,  as  shown  in  Figure  2.5 
(b).  Assuming  (for  the  moment)  that  no  poles  or  branch  points  lie 
within  the  enclosed  region,  the  following  relations  hold  from 
fundamental  principles: 


(2.81) 


where  C-  is  a  path  along  the  real  axis;  C  is  a  circular  arc  of 
i  R 

arbitrarily  large  radius  R;  and  is  the  actual  transformed  contour 
(from  the  original  real  £-axis).  So  long  as  the  integral  over  CR 
vanishes,  the  original  real-axis  integration  in  terms  of  £  ,  Equation 
(2.57),  can  again  be  expressed  as  a  real-axis  Integration  in  terms  of 
t.  That  is, 


ive;  and 
Quadrant 
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£“°°  £  MOO 

|  dC  -  |  dt  -  |  dt  -  |  dt  .  (2.82) 

C"o  Ct  C x  t-o 

As  this  is  the  desired  result,  the  final  step  in  the  transformation  is 
to  ensure  that  the  integral  over  CR  vanishes  as  R  gets  large. 

The  exponent  in  Equation  (2.76)  can  be  written  in  terms  of  its 
real  and  imaginary  parts  for  and  B^B^+iB^  as18 

-ikR2(t2  +  2Bt) --ikR2[(tr  +  iti)2  +  2(Br  +  iBi)(tr  +  iti)]  ,  (2.83) 

where  the  real  part 

Re[-ikR2(t2  +  2Bt)]  -  2kR2trti  +  2kR2Brt±  +  2kR2Bitr  (2.84) 

must  be  less  than  zero.  First,  for  path  CR  "starting"  on  the  real 
axis,  t j»0,  so  that  ReI-ikR2(t2+2Bt)]  -  2kR2Bitr.  Hence,  to  ensure 
convergence  here,  B^  must  be  negative.  Furthermore,  for  all  other 
points  along  CR  off  the  real  axis,  the  first  term  on  the  right  hand 
side  of  Equation  (2.84)  is  sufficiently  large  to  ensure  that  the 
exponent  is  negative  regardless  of  the  sign  of  Br.  Consequently,  the 
proper  root  for  B  is  now  fixed,  and  Equation  (2.65)  may  be  re-written 
to  enforce. the  constraint  B^<0  as 

18 

The  wavenumber  k,  which  had  been  assumed  complex  for  physical 
purposes,  can  here  be  assumed  real  for  mathematical  purposes.  In  the 
alternative,  the  small  positive  Imaginary  part  k  could  be  made 
arbitrarily  small  so  that  k  ■  kf  +  ikt  — ►  kp  in  the  limit. 


-i  A  + 


$slni|> 


$  cost|> 


Re  /"  >0 


(2.85) 


and  Che  transformed  path.  Equation  (2.77)  can  be  specified  as 


t(€) 


>  0 


(2.86) 


where  the  symbol  indicates  that  the  root  is  to  be  taken  to  yield  a 
positive  imaginary  part. 


2.3.3  The  Branch  Points  and  Branch  Cuts 

The  simple  pole  that  was  present  in  the  original  integral  (see 
Equation  [2.48])  does  not  appear  in  the  transformed  integral  of 
Equation  (2.76);  however,  due  to  the  multivalued  nature  of  the 
denominator,  a  branch  point  is  still  present.  It  can  be  shown  that 
the  presence  of  a  branch  cut  arising  from  this  point  will  not  affect 
the  integration  directly  since  the  contour  can  be  closed  without 
crossing  any  portion  of  the  cut.  But,  as  will  be  pointed  out  later, 
the  branch  point  is  also  a  singularity  of  the  integrand  and,  as  such, 
will  limit  the  radius  of  convergence  of  the  Taylor  series  expansion  of 
the  denominator. 

The  branch  points  are  simply  the  roots  of  the  equation 


t 

"h 


2Bt 

H 


-  -  -sat  +  ! 


(2.87) 


where,  again, 


H  ■  1  +  BsiniJj  +  /I  -  p  co sip 


(2.88) 


Therefore,  the  branch  points  are 


where 


t^  ■  -B  ±  /2A*  ;  Re  >  0 
2 

A  ■  A  -  B2'  cosi^  . 


(2.89) 


The  branch  cuts  associated  with  these  points  have  already  been  chosen, 
implicitly.  That  is,  the  specification  that  the  root  in  the  integral 
of  Equation  (2.76)  be  taken  to  yield  a  positive  real  part  has 
essentially  defined  lines  in  the  t-plane,  across  which  the  value  of 
the  root  changes  in  sign.  The  equation  describing  these  branch  lines 
can  be  derived  in  a  relatively  straightforward  manner. 

In  the  complex  z-plane,  the  negative  real  axis  acts  as  the  branch 
cut  for  the  constraint  Re(z)**>0;  that  is,  the  equation  describing  the 
branch  line  is  z«.  ■  -p,  for  p£0.  In  direct  analogy,  the  constraint 


>  0 


(2.90) 


translates  into  the  defining  equation  for  the  branch  line 
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Solving  this  for  t,  gives  the  equation  for  the  branch  line  in  the  t- 
plane.  Thus, 

tllne  -  -B  ±  /2A  +  pH'  ;  p*0  .  (2.92) 

It  can  be  shown  that  the  branch  lines  described  by  this  equation 
usually  lie  wholly  within  the  first,  second,  or  third  quadrants  of  the 
t-plane,  and  consequently  do  not  interfere  with  the  closed  integration 
contour.  For  the  rather  unrealistic  case  for  which  B^O,  |B|  »1,  and 
slnij/^O,  one  of  the  branch  points  can  fall  in  the  fourth  quadrant; 
however,  the  corresponding  cut  extends  downward,  and  the  closed 
contour  can  be  perturbed —  if  necessary— >to  avoid  crossing  it. 

The  transformed  path,  Equation  (2.86),  along  with  the  branch 
cuts.  Equation  (2.92),  are  plotted  in  Figure  2.6  (a)  and  (b)  for  two 
representative  combinations  of  admittance  B  and  grazing  angle  The 
contours  are  presented  here  only  to  add  insight  into  the  problem,  as 
the  final  integration  need  only  be  performed  over  the  positive  real 
axis. 

2.3.4  The  Taylor  Series  Expansion 

The  Integral  in  Equation  (2.76)  can  be  re-written  now  as  a  real 
Integral.  Thus,  disregarding  the  constant  multiplier, 
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I 


3 


4 


-ikR2[t2  +  2Bt] 


e 


dt 


(2.93) 


The  square  root  function  in  the  integrand  is  still  the  obstacle 
preventing  the  evaluation  of  this  integral,  so  the  radical  will  be 
removed  by  expanding  it  into  a  Taylor  series  about  the  origin.  It  is 
possible  to  avoid  computing  the  derivatives  normally  required  by 

Taylor's  formula  by  noting  that  the  given  function  is  already  in  the 

2 

standard  form  (1+u)  where  u*-(t  /H  +  2Bt/H).  From  the  well-known 
result 


(1  +  a)-1*  -  1  -  -J  »  +  |  u2 


5  3  35  4 

16  u  128  u 


(2.94) 


the  expansion  for  the  integral  factor  can  be  translated  to  be 
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(2.95) 


n“0 


A  general  form  for  the  Taylor  coefficients  Tr  can  also  be  obtained  as 


follows: 


(2.96) 


where  the  quantity  j  is  the  binomial  coefficient;  the  parameter  G 

is  given  by 


1  +  Ssinip  - 


costjj  ; 


(2.97) 


and  a  ^  is  expressed  by  the  following  recursion  relation: 


(2.98) 


The  notation  on  the  summation  symbol  in  Equation  (2.96)  indicates  that 

k  increases  in  Integer  steps  until  the  expression  (n-2k)>0  is  no 

longer  satisfied  (for  example,  if  n«6  or  n-7,  then  k“0,l,2,3).  For 

reference,  the  values  of  the  binomial  coefficient  and  a  required  for 

01 

an  expansion  through  t^  are  given  as  follows: 


afi  -  231/1034 


-  20 
-  15 


As  with  any  Taylor  series,  Che  expansion  in  Equation  (2.95)  is 
absolutely  and  uniformly  convergent  within  any  circular  region  in  the 
complex  plane  about  the  expansion  point  (t"0  in  this  case)  that  does 
not  contain  a  singularity.  The  circle  of  convergence  for  Equation 
(2.95)  is  thus  limited  by  one  of  the  branch  points  t^  or  defined  in 
Equation  (2.89).  Precisely,  if  |ta|  represents  the  smaller  of  the 
magnitudes  |t^|  and  Itjl*  then  the  radius  of  convergence  for  the 
Taylor  series  will  be  |t  |.  Of  course,  for  values  of  t  such  thar 
|t|>|tal,  the  function  on  the  left  hand  side  of  Equation  (2.95)  can  no 
longer  be  equated  to  the  series  on  the  right,  because  the  latter  may 


diverge 


66 


2.3.5  The  Term-by-Term  Integration 


The  integral  in  Equation  (2.93)  can  now  be  re-written  as  a  sum 
of  two  convergent  integrals 


t'  +  T" 

i3  L3 


p  -ikR2(t2  +  2Bt)  ” 


Zvn 


®  -ikR2(t  +2Bt) 
dt+  f  - - dt  . 


Ital  A  _  tl.  2Bt 

H  H 


(2.100) 


It  should  be  emphasized,  here,  that  no  mathematical  approximations 
have  been  made  thus  far  in  the  derivation  and  so  the  exact  expression 
for  the  field.  Equation  (2.47),  can  be  written 


ikR  ikR.  ikR, 

..  .  e _ _  ,  e _ ‘  4ikRBBe  . 

*(r,Z/  R  R  "  (B  +  l3  ' 


(2.101) 


with  I ^  given  above. 

To  insure  accuracy  in  the  approximations  to  be  made,  the  relative 
magnitudes  of  the  two  Integral  terms  in  Equation  (2.100)  should  be 
investigated.  If  the  second  integral  can  be  neglected,  then  a  term- 
by-term  integration  on  the  first  integral  would  accurately  represent 


the  field.  It  is  sufficient  to  examine  first  the  behavior  of  the 


Integrand 
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-ikR2(t2  +  2Bt) 

U(t)  =  - -  .  (2.102) 


Plotted  In  Figure  2.7  (a),  (b),  and  (c)  are  the  real  and  imaginary 

parts  of  the  function  U(t)  versus  t  for  three  values  of  the  parameter 

kR2«  Here,  3  and  ip  are  fixed  at  the  representative  values 

0“(O .055 ,-0 .062)  and  ^3.8°.  The  value  of  |t  |  for  this  B  and  is 

a 

marked  by  a  vertical  bar  on  the  t-axis  and  indicates  the  limit  for  t 

beyond  which  the  Taylor  series  diverges*  It  is  apparent  from 

Figure  2.7  that  as  kR2  increases  a  greater  and  greater  contribution  to 

the  integral  comes  from  the  region  in  which  the  Taylor  series 

converges.  In  other  words,  for  large  kR2>  the  integral  from  t*0  to 

t«°°  in  Equation  (2.93)  can  be  replaced  by  just  the  first  integral  in 

Equation  (2.100).  The  condition  "large  kR2"  can  often  be  taken  in  a 

relaxed  sense.  In  the  particular  case  shown  if  Figure  2.7  (b),  for 

Instance,  the  relatively  modest  value  of  kR2**13.5  shows  a  very  minor 

contribution  for  t^lt  I. 

a  1 

There  Is  one  final  step  to  be  taken  before  the  integral  in 
Equation  (2.93)  can  be  expressed  in  a  form  suitable  for  term-by-term 
integration.  This  involves  replacing  the  denominator  in  the  second 
integral  of  Equation  (2.100)  by  the  Taylor  series  expansion,  and  thus 
writing 
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I 


4 


oo  oo  9 

£  T  J  t'  +  210 


(2.104) 


Clearly,  from  the  discussion  in  the  previous  section,  the  second 

Integral  in  Equation  (2.103)  does  not  equal  the  second  integral  in 

Equation  (2.100)  and  so  1^  j4  1^.  However,  a  relationship  does  exist 

between  the  two  integrals,  namely  that  1^  given  by  Equation  (2.104)  is 

an  asymptotic  expansion  (sometimes  called  a  ”semi-convergent" 

expansion)  of  the  integral  I^.  An  asymptotic  expansion  converges  up 

to  a  point,  and  then  diverges.  That  1^  represents  an  asymptotic 

expansion  of  1^  is  proven  in  principle  in  several  classical 
19 

mathematics  texts  and  may  be  explained  heuristically  here  as 
follows . 

The  second  series  of  Integrals  on  the  right  hand  side  of  Equation 
(2.103)  clearly  diverges  since  it  involves  the  Taylor  series  outside 
of  its  circle  of  convergence,  but  the  first  series  converges.  The 


1 Q 

The  necessary  tools  are  Watson's  lemma  and  some  theorems  on 
Laplace  transforms.  See,  for  example,  Copson,  1935,  p.  218. 


rates  of  convergence  and  divergence  are  controlled  by  the  parameter 
kR2  in  the  exponent:  the  larger  kR2,  the  more  rapid  the  convergence 
of  the  first  series,  and  the  slower  the  divergence  of  the  second. 
However,  for  any  fixed  value  of  kR2,  the  sum  of  the  two  integral 
series  behaves  asymptotically  since  it  will  always  converge  up  to  a 
certain  number  of  terms,  at  which  point  the  second  series  "takes  over" 
and  causes  it  to  diverge.  Moreover,  recalling  the  conclusions  drawn 
from  Figure  2.7  (a),  (b),  and  (c),  the  series  up  to  the  onset  of 
divergence  is  an  accurate  representation  of  the  exact  integral  for 
"large  kR2.“ 

The  conclusion  to  be  drawn  from  this  section  is  that 


I3- 


•  -ikR2(t  +  2Bt) 
e _ 


* 

4 


dt 


j  / - j ' 

o  J  t  2Bt 


-Z  T»l  ** 


oo  o 

-ikR  (t  +  2Bt) 

e  ^  dt  .  (2.105) 


n*o 


H  H 


where  the  symbol  translates  "has  the  asymptotic  expansion"  or  "is 
asymptotically  equal  to."  Each  Integral  in  the  series  is  tractable, 
and,  as  a  result,  it  will  be  possible  to  express  the  total  field 
solution  as  an  asymptotic  series  of  closed-form  terms. 


2. 3. 5.1  The  Formal  Series  Solution.  The  general  integral  in  the 


series  Equation  (2.105)  can  be  solved  in  terms  of  parabolic  cylinder 
functions.  The  functional  form  is  (Gradshteyn  &  Ryzhik,  1965, 
Integral  #3.462[1]): 
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n  -at  -bt 
t  e  at 


(2a) 


IiH-iy2r<°+l)  eb2/8‘  D-(.*1)(b"K> 


(2.106) 


where  r(m)  Is  the  gamma  function— for  integer  arguments  r(m-fl)»m! 

— and  D_^(z)  is  the  parabolic  cylinder  function.  Several  fundamental 
relationships  involving  the  functions  D_^(z)  will  be  useful  for  later 
development  (see  Gradshteyn  &  Ryzhik,  1965,  pp.  1064-1068).  First, 


D_^(z)  -  t/j  eZ  ^  erfc(z//2  ) 


(2.107) 


and 


20 


D_2(z) 


ez  ^  erf  c  (zh/l  )  -  e  2  ^  , 


(2.108) 


where  erfc(z)  is  the  complex  complementary  error  function  defined  as 


erfc(z) 


w 

X  [  e-T: 


dT  . 


(2.109) 


In  addition,  parabolic  cylinder  functions  obey  the  following  recursion 
relationships: 


20 

There  is  a  sign  error  in  the  expression  for  D_  (z)  found  in 
the  text  by  Gradshteyn  and  Ryzhik.  The  expression  givln  here  is  the 
corrected  form. 
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or,  in  general, 
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rtr^r^W  ■■  J»".»  rr  '.9  TT*  _■  g*  .■  .»  .■  ■  .■  »-  ~  .  r.  f  ;T;  ,  i-w-  -  i  -  *—  . 
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Noting  that  r(m+l)-(m+l)r(m) ,  the  coefficients  E  and  K  can  be 

ID  01 

expressed  more  concisely  as 


E 

m 


_b_ 

2a 


and 


_b_ 

2a 


(2.120) 


(2.121) 


where 


E  -  1  K  "  0 

o 

(2.122) 

-  b  -  l 

1  2a  K1  2a  * 

The  general  integral  at  the  beginning  of  this  section  has  thus 
been  solved  in  terms  of  the  complementary  error  function  and  two 
general  coefficients,  for  which  recursion  relations  are  provided.  The 
specific  integral  of  interest  here— Equation  (2.105) — can  now  be 
expressed  from  the  above  results  by  translating  ”a"  to  "ikR,^  and  "b" 
to  2ikR2B.  Thus,  the  asymptotic  expansion  of  the  integral  in 
Equations  (2.93)  and  (2.105)  can  be  written  in  the  final  form 

00 

X3  -  E  Tn  [6o  En  +  Kn]  (2.123) 

n«o 


where 
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E 

m 


iszil 

2ikR2 


•  -“.-i +  k.-2 


(2.124) 


All  complex  square  roots  in  the  above  equations--including 
2  k 

(1-6  are  to  be  taken  to  yield  a  positive  real  part.  It  must  be 
stressed  that  the  definition  for  X  in  Equation  (2.127)  has  been 
derived  from  the  argument  of  the  parabolic  cylinder  function,  and  the 
roots  must  be  taken  separately  as  Indicated.  This  is  Important 
because,  in  general,  /ipf  (z)  +  )/Tp  >/f  (z) ,  when  p  is  a  positive 
constant  and  positive  real  roots  are  taken. 
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The  full  solution  for  the  field  at  the  point  (r,z)  can  now  be 
written  as 


t(r,z) 


IkR  ikR2 


ikR, 


4ik8Be 


(6  +  sini^) 


I 


n»o 


T  [e  E  +  K  ] 
non  n 


(2.128) 


I 

I 
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2. 3. 5. 2  The  Q-term  Form  of  Solution.  As  mentioned  in  section 
2.1,  it  is  often  desirable  to  express  the  total  field  as  a  simple  sum 
of  a  direct  component  from  the  source  plus  a  reflected  component  that 
appears  to  originate  from  an  "image”  source  at  (r,-z^).  Hence,  a 
solution  of  the  form 


<Kr,z) 


+  Q 


(2.129) 


will  be  sought  in  this  section. 

The  last  term  in  Equation  (2.128)  can  be  re-written  as 


6  e 


ikR.  £  08 


3  (6  +  sin40  R2 


Z[-2ikBR„e  T  E  -2ikBR  T  K  ]  .  (2.130) 

2  o  n  n  2  n  n 


n-0 


Furthermore,  it  is  convenient  to  re-express  the  product  TnEn  as 
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(2.134) 
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When  Equation  (2.130)  is  substituted  back  into  Equation  (2.128), 
the  total  field  may  be  expressed  in  the  form  of  Equation  (2.129),  with 
the  "spherical  wave  reflection  coefficient"  Q  given  by 


Q 


1  + 


28 


(8  +  siniJO 


I 


n=o 


*  *  * 

T  [e.E  +  K  ] 
n  1  n  nJ 


(2.135) 


Here,  the  parameter  e^  has  been  introduced  where 


e  =  -2ikBR0e  ■  -kR.  \/G  /t?'  e  ^  erfc(-iX) 

1  L  O  L  lKK.^ 


i/ir  X  e  ^  erfc(-iX)  , 


(2.136) 


and  X  is  defined  by  Equation  (2.127). 

2. 3. 5. 3  The  F-term  Form  of  Solution.  The  asymptotic  series 
solution  developed  in  section  2.3.5. 1  can  also  be  put  into  a  form  that 
provides  a  more  direct  comparison  to  the  classical  plane  wave 
solution.  As  such,  the  solution  would  comprise  a  direct  component 
from  the  source,  a  "plane-wave  reflected"  component,  plus  a  correction 
term  that  accounts  for  the  sphericity  of  the  wavefronts  (see  Equation 
[2.7]).  Historically,  several  approximate  solutions  in  both 
electromagnetics  and  acoustics  have  been  expressed  in  such  a  form,  and 
these  contain  the  so-called  ’F-term"  or  "attenuation  factor."  The 
solutions  are  expressed  as 
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ikR  ikR2  ikR2 

$(r,z)  =  - -  +  R  -^-r—  +  (1  -  R  )  F  ^-r— 

R,  P  R-,  P  R-, 


(2.137) 


where  is  the  plane  wave  reflection  coefficient. 


__  sinijj  -  8 

p  sinij;  +  8 


v  2 . 1 38  ) 


This  form  of  solution  is  easily  derived  from  the  Q-term  form 
given  in  the  previous  section.  Thus, 


28 


(8  +  sin({0 


sinij;  +  8  -  sinijj  +  3  =  . 

(8  +  sinty)  P 


and,  therefore 


Q=*Rp  +  (l~Rp)  F 


E*  *  * 

T  [e  E  *  K 
n  1  n  r. 


n»o 


Solving  this  for  F  gives  the  final  rrs 


F  «  1  - 


I  <  ' 
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2. 3. 5. 4  The  Hankel  Function  Form  of  Solution*  The  results  of  at 
least  one  recent  study  (Chien  &  Soroka,  1975),  which  applied 
asymptotic  techniques  to  the  point  source  propagation  problem,  were 
expressed  directly  in  terms  of  the  zero-order  Hankel  function.  The 
same  function  has  also  appeared  in  component  terms  (most  often,  the 
"surface  wave"  term)  in  other  solutions  offered  in  the  literature.  In 
view  of  this,  the  asymptotic  series  solution  will  be  expressed  in  an 
alternate  form  that  is  also  in  terms  of  the  Hankel  function.  It  is 
offered  here  mainly  to  facilitate  comparisons  with  other  solutions,  as 
the  series  forms  presented  so  far  in  this  study  are  the  preferred 
forms  of  solution.  (The  series  solutions,  in  conjunction  with  the 
recursion  formulae,  translate  readily  into  computer  programs  for 
providing  numerical  data.) 

The  first  few  product  terms  T*E*  from  the  Q-term  solution. 
Equation  (2.135),  can  be  written  out  explicitly  to  generate  the 
following  series: 


83 


where  Che  parameters  H  and  G  are  repeated  here  as 


H  -  1  +  BsiniJ/  +  (1  -  V  costjJ 
G  -  1  +  Bsini|>  -  (1  -  02)1*  cosip 


(2.143) 


Each  of  the  series  In  the  brackets  subscripted  a,  b,  and  c  can  be 
recognized  to  be  simple  expansions  of  the  form  (l+x)a  such  that 

r  r  ^-1/2 


[]. 

U 

[1 


MI- 

MI 

MI 


(2.144) 


Using  the  fact  that  [  1— (G/H) J  »  [(F-G)/HJ  -  2A/H,  where  A-O-g^cos^, 
the  partial  series  in  Equation  (2.142)  can  be  re-written  as 


Z  TX  ■  fi 


Partial 


fjL 

fi  +  1  f- 

F 1  ,  9  Tf| 

W 

A  4ikR2F^ 

^  32(ikR2)F2'-2AJ 

(jl 

H 

r,  .  i  . 

9  i  1 

[lAj 

1  8ikR2A 

2  + 

128(ikR2A)^ 

Comparing  this  to  the  asymptotic  expansion  of  the  Hankel  function 
H^^(z)  (see,  for  example,  Abramowltz  &  Stegun,  1964,  Section  9.2.8) 


(2.146) 


bo1><2>  •  •“  [i + is + 


128 (iz)' 


-] 


enables  Equation  (2.145)  to  be  re-phrased  as 


i « -  W 


-ikR^A  q. 

e  L  Hy;(kR,A)  . 
o  i. 


(2.147) 


n*o 


With  the  added  fact  that  HG"(S+slniJ>)  ,  the  total  solution  given  by 
Equations  (2.129)  and  (2.135)  can  be  expressed  In  the  final  form 


ikR  lkR 

„  1  m  2 

Hr,z)  -  -S-g—  +  —— 

R1  R2 


-lkR  Bsinl|»  ... 

-  kfiire  z  erfc(-iX)  Hu;[kR,(l 

O  L 


B2)1*  cosi|»] 


(2.148) 
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2B  e _ _  *  * 

+  (B  +  slntfO  R.  L-t  ^nRn 


n-o 


2. 3. 5. 5  The  Solution  Using  the  First  Term  Only.  For  many  cases 

2i 

of  practical  Interest,  the  parameter  kR^  is  large  and  the  asymptotic 
solution  Is  rapidly  convergent  In  the  initial  terms.  Thus,  using  only 


21 

As  examples:  At  1000  Bz.,  for  source  and  receiver  on  the 
ground  and  separated  by  100  feet,  the  value  of  kR.  is  559.  At  200  Hz, 
for  source  and  receiver  4  feet  above  the  ground  and  separated  by  20 
feet,  the  value  of  kRj  Is  22. 


the  first  (n-0)  tens  in  Equation  (2.135)  often  provides  sufficient 
accuracy.  In  such  instances,  the  total  solution  can  be  written  as 
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ikRL  ikR2 

where 

-X2 

■  1  +  i/t  X  e  erfc(-iX) 


(2.149) 


(2.150) 


(2.151) 


The  Solution  for  Perpendicular  Incidence 


Although  the  asymptotic  series  solution  is  valid  for  ^*90° 
(receiver  directly  above  or  below  the  source),  an  exact  solution  can 
be  obtained  for  this  special  case  in  a  straightforward  way.  Thus, 
when  slnij;  -  1  is  substituted  into  the  original  transformed  integral. 
Equation  (2.67),  and  it  is  noted  that  (1  -  £2  +  2i£)  ■  (1  +  i£)^,  the 
integral  can  be  re-written  as 
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•  ikR2[l  +  i(l  +  6)51 


J-^rr 


21k8e 


11R,  7  -»*2a  +  «« 

2  j  i~a- 1)  « 


(2.152) 


Using  the  relation  (Gradshteyn  &  Ryzhlk,  1965,  Integral  #3.352[4J) 


J  dx  -  -  eBuEi(-yS) 


[|argB|<ir  ;  Rey>0] 


(2.153) 


where  El  la  the  exponential  Integral  defined  by 


El  (w) 


w 

(  x 

I  irix  > 


(2.154) 


the  Integral  1^  In  Equation  (2.152)  becomes 


-ikILB 

Ix  -  -21kBe  EltikR2(l  +  3)]  . 


(2.155) 


Furthermore,  El(w)  has  the  convergent  series  expansion  (Abramowltz  & 
Stegun,  1964,  Section  5.1.10) 


Ei(-w)  - 


y  +  lnw  +  ^ 
n“l 


(-1)°  wn 


(2.156) 


where  Y  is  the  Euler  constant,  Y"0. 57721  56649...  The  exponential 
Integral  also  has  the  asymptotic  expansion  for  large  argument 
(Abramowltz  &  Stegun,  1964,  Section  5.1.51) 
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W  1  W  w2  w3  -/ 
00 

e”W  V  (-1)^  k!  .  i  i  .  3ir 

—  )  * — ^ J  |arg  w|  <  — 


(2.157) 


In  summary,  the  total  field  at  the  receiver  point  for 
perpendicular  incidence  (r-0)  Is 


ikRL  ikR2  .ikR  g 

•  (*)-■——  +  ■— - 2ik8e  2  Ei[ikR,(l  +  B)]  ,  (2.158) 


in  which  R2  nay  be  replaced  by  R2  »  h  »  z  +  zQ. 

When  kR2  is  very  large,  the  first  term  in  the  asymptotic  series. 
Equation  (2.157),  may  be  used  so  that  the  total  field  becomes 


l2  ”  R2(l  +  B) 


ikR, 

±_i+  fi-  e| 

*!  U  +  ej 


(2.159) 


thus  reducing  to  the  "plane  wave"  solution  (Equations  [2.4]  end  [2.5] 
with  sin  \f>  ■  1). 


CHAPTER  III 


A  MODEL  FOR  THE  COMBINED  PROBLEM  OF 
POINT  SOURCE  PROPAGATION  WITH  A  DIFFRACTING  BARRIER 

When  a  barrier  located  on  an  Impedance  plane  obstructs  the  line 
of  sight  from  a  point  source  to  a  point  receiver,  acoustic  energy 
reaches  the  receiver  by  diffraction  over  the  edge  of  the  barrier.  The 
"canonical"  problem  of  point  source  diffraction  by  an  ideal  half  plane 
is  relatively  well  understood;  however,  the  attempts  to  apply  its 
solution  directly  to  the  problem  here  have  met  with  little  success. 
This  is  due  to  the  fact  that  when  a  barrier  is  located  on  the  ground, 
the  incident  field  at  its  edge  is  not  a  simple  monopole  field. 

Instead,  it  is  a  complex  field  resulting  from  the  propagation  over  the 
Impedance  plane.  Furthermore,  the  ground  plane  interactions  on  the 
receiver  side  of  the  barrier  must  be  accounted  for  when  predicting  the 
total  field  behind  the  barrier* 

In  this  chapter,  an  "Edge-Plus-Images"  model  will  be  described 
that  merges  the  propagation  results  of  the  preceding  chapter  with 
half-plane  diffraction  theory.  The  mathematical  model  is  readily 
programmable  on  a  digital  computer,  and  thus  practical,  engineering 
calculations  can  be  obtained  with  little  trouble.  Although  the 
following  discussion  is  in  terms  of  a  vertical  barrier  on  a  ground 
plane,  the  results  are  general  and  may  be  applied  to  any  thin,  planar 
protrusion  from  a  flat,  impedance-covered  surface. 


89 


3.1  The  Geometrical  Theory  of  Diffractloon 


The  use  of  Keller's  Geometrical  Theory  of  Diffraction  (GTD)  in 
solving  practical  problems  has  met  with  universal  success,  and  it  will 
serve  here  as  the  foundation  in  building  the  diffraction  model  for  a 
barrier  on  a  ground  plane.  The  derivation  of  the  GTD  formulation  will 
not  be  given  here  since  this  information  is  available  in  many  texts  as 
well  as  in  Keller's  original  papers  on  the  subject  (see  the  references 
cited  in  Chapter  1).  The  predictions  of  the  theory  have  been  tested 
and  confirmed  many  times  in  both  the  fields  of  electromagnetics  and 
acoustics.  Agreement  with  exact  analytical  solutions,  when  available, 
has  been  excellent;  correlation  with  experiments,  in  the  absence  of 
mathematical  solutions,  has  also  been  encouraging.  Furthermore, 
although  it  is  in  principle  a  high-frequency  method,  the  GTD  often 
performs  with  a  high  degree  of  accuracy  down  to  wavelengths  comparable 
to  the  size  of  the  scattering  object. 

Essentially,  the  GTD  supplements  the  approach  of  dassslcal 
geometrical  acoustics  by  including  "diffracted  rays,”  which,  as 
opposed  to  "specularly  reflected  rays,"  can  penetrate  into  the  "shadow 
zone."  And  just  as  the  law  of  conservation  of  energy  flux  in  a 
"bundle  of  rays"  (Fermat's  principle)  leads  to  Snell's  law  for 
specular  reflection  (angle  of  incidence  »  angle  of  reflection), 
similar  arguments  applied  to  the  diffracted  rays  lead  to  the  "law  of 
edge  diffraction."  With  reference  to  Figure  3.1,  Fermat's  principle 
prescribes  that  the  total  path  length  for  the  edge-diffracted  ray  S-E- 
P  should  be  a  minimum  (out  of  all  possible  paths  connecting  the 


Figure  3.1  The  Law  of  Edge  Diffraction:  Incldeni 
Diffracted  Ray  Subtend  Equal  Angles  0 
Edge.  Diffracted  Rays  form  a  Cone  wii 
Diffraction  Point  E. 
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points)*  Hence,  the  law  of  edge  diffraction  can  be  stated  as:  when 

an  incident  ray  strikes  an  edge  at  an  angle  0  ,  the  set  of  diffracted 

rays  will  comprise  all  those  rays  that  also  make  an  angle  0  (in  the 

22 

"forward’*  direction)  with  the  edge*  The  resulting  "cone"  of 

diffracted  rays  is  illustrated  in  the  figure*  If  the  edge  is  straight 

(as  it  is  in  the  present  barrier  problem),  the  point  E  will  be  unique, 

23 

and  only  a  single  diffracted  ray  will  reach  the  receiver. 

Generally,  the  edge-diffracted  field  is  described  in  the  GTD 
formulation  by  several  factors:  (a)  a  reference  field  at  some  point 
along  the  ray,  (b)  a  geometrical  spreading  factor  depending  on  the 
radii  of  curvature  of  the  edge  and  of  the  incident  and  diffracted  ray 
bundles,  and  (c)  a  phase  factor  accounting  for  the  propagation  of 
energy  over  the  entire  length  of  the  ray  path.  The  reference  field  is 
usually  given  at  the  point  E  on  the  edge  and  is  expressed  in  terms  of 
a  "diffraction  coefficient,"  which  describes  the  ratio  of  the 
diffracted  energy  to  the  incident  energy.  It  can  be  shown  that  for  a 
spherical  wave  incident  upon  a  straight  edge,  the  diffracted  field  at 
the  receiver  point  F  can  be  expressed  in  GTD  terms  as: 


22 

The  notation  ”0”  for  this  angle  has  been  adopted  here  in  order 
to  be  consistent  with  the  customary  usage  in  the  literature. 

Hopefully,  the  context  will  prevent  any  confusion  with  the  admittance 

"0." 

23 

If  the  edge  were  not  straight  as  shown,  but  curved,  the  angle 
0  would  be  measured  from  the  line  tangent  to  the  edge  at  the  point  E. 
Thus,  the  possibility  for  several  rays  reaching  the  receiver  from 
different  points  along  the  edge  would  exist. 
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♦dlff(?)  '  *inc 


/ - R - V 

«>  ».<♦.*♦•»  4,(8,  '»2) 


ikR„ 


(3.1) 


where  the  angles  and  distances  are  defined  as  in  Figure  3.2.  Here, 
$inc(E)  is  the  incident  field  at  the  point  E  on  the  edge.  For  a 
simple  point  source,  therefore. 


ikR, 


*i»<®  • 


(3.2) 


The  factor  Dg  is  the  diffraction  coefficient  for  a  straight  edge;  the 
square  root  quantity  is  the  geometrical  spreading  factor;  and  the  last 
term  accounts  for  the  phase  along  the  ray  from  the  edge  to  the 
receiver. 

Probably  the  strongest  advantage  of  the  GTD  formulation  is  that 
all  the  relevant  properties  of  the  diffracting  obstacle  are  contained 

in  the  diffraction  coefficient  itself.  If  a  wedge  were  being 

✓ 

considered  instead  of  a  straight  edge,  the  diffracted  field  would  have 
the  same  form  as  Equation  (3.1)  only  the  diffraction  coefficient  would 
be  derived  from  wedge-diffraction  theory.  For  objects  with  complex 
geometries,  the  diffraction  coefficient  would  be  taken  from  the 
canonical  problem  (half-plane,  wedge,  sphere,  cylinder,  etc.)  that 
most  closely  describes  the  local  shape  of  the  object  at  the 
diffraction  point. 
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In  Keller's  original  work  on  the  Geometrical  Theory  of 
Diffraction  applied  to  the  rigid  or  pressure  release  half-plane 
problem,  he  deduced  the  diffraction  coefficient  from  the  first  term  in 
the  expansion  of  Sommerf eld's  exact  solution.  Thus, 


-,  rigid 


+,  pressure  release 

The  diffraction  coefficient  has  singularities  at  <f>  «  -tt+  <f>Q  and  at 
$  -  ir  -  <J>q,  which  correspond  to  the  physical  shadow  boundaries  as 
shown  in  Figure  3.3.  Furthermore,  when  sing  ->0,  the  solution  becomes 
invalid.  This  latter  condition  corresponds  to  zQ  in  Figure  3.3  being 
very  large  relative  to  r^  and  r. 

Many  researchers  have  modified  the  Sommerfeld  coefficient  or 
derived  Independent  expressions  in  terms  of  Fresnel  Integrals,  which 
remain  valid  in  the  transition  regions  of  the  shadow  boundaries.  The 
more  exact  diffraction  coefficient  for  a  rigid  half-plane  can  be 
written  (Hayek,  in  Varadan  &  Varadan,  1982)  as 


Figure  3.3  Angles  Defining  Incident  and  Reflection  Shadow 
Boundaries . 
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where  sgn(x)  means  the  sign  of  the  argument  x,  and 
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and  F* (x)  -  C(x)  +  i  S(x)  , 


where  C(x) 
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cost  dt  ,  S(x) 
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are  the  Fresnel  cosine  and  sine  Integrals,  respectively. 


3.1.2  The  GTD  Solution  for  an  Inpedance  Covered  Half-Plane 


Recently,  an  exact  Integral  solution  has  been  derived  (Kendig, 
1977)  for  the  diffraction  by  a  half-plane  with  locally-reacting 
impedance  surfaces.  The  advantage  of  this  solution  over  a  similar 
contribution  by  Halyushinets  (1955,  1962)  is  that  Kendlg's  exact 
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Integrand  is  in  closed  form,  thus  allowing  asymptotic  techniques  to  be 
applied  to  the  integral  for  its  evaluation.  The  resulting  solution 
has  been  adapted  to  the  GTD  formulation  (Hayek,  et  al.,  1978)  and  the 
diffraction  coefficient  can  be  written: 


, _  -ikA 

sgn  (A)  4.(4)  JkL  e  _±/u 

D  (<>,*  ,6)  * - TiTH -  [1  "  &  e  F* ( | A | ) ] 


e  o 


sin  8 


sgn  (B)  4j(4)  /kL  e 
sin  8 


(3.6) 


[1  -  ^2  e“i7r/4  F*(|b|)] 


i 

in  which  the  parameters  are  defined  by  Equation  (3.5).  In  addition, 
the  special  functions  4^  and  4^  contain  the  Impedance  Information  and 
can  be  expressed  as 
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The  impedance  of  Che  half-plane  surfaces  has  been  expressed  in  terms 
of  the  so-called  "Brewster"  angles  0*" ,  for  the  face  at  ♦  -  +7,  and  0”, 
for  the  face  at  $  ■  -7,  where 


sine*  -  -r22 - 

Z  sing 


(3.9) 


Because  the  field  at  the  receiver  is  affected  more  by  the 
Impedance  on  the  lnsonlfled  (source-facing)  surface  of  the  half-plane 
than  by  that  on  the  "dark"  surface,  an  expression  for  the  case  when 
0~  **  0  (rigid)  and  6+  arbitrary  might  be  useful  for  practical  cases. 
Thus,  for  the  impedance-rigid  half-plane,  the  special  functions 
become : 
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3.2  The  GTD  Model  for  the  Barrier  on  a  Ground  Plane 

Having  defined  the  edge— diffraction  coefficients  to  be  employed 
in  the  Geometrical  Theory  of  Diffraction  formulation,  it  is  now  a 
simple  matter  of  incorporating  the  ground  effects  into  the  model. 

Just  as  the  diffraction  effects  were  lumped  into  the  coefficient 
D(*0,<M)  in  Equation  (3.1),  the  incident  field  at  the  barrier  edge  is 
totally  contained  in  the  term  4^nc(E) .  When  a  locally-reacting  ground 
plane  is  present  beneath  the  source,  *inc(E)  is  no  longer  a  simple 
monopole  field  but  instead  represents  the  total  field  that  exists  at 
the  edge  point  E.  This  field  is  derived  from  the  ground  propagation 
solution  described  in  Chapter  2  and  will  be  modelled  as  a  super¬ 
position  of  two  ray  paths — the  direct  path  and  the  ground-reflected 
"image*'  path.  Similarly,  the  energy  spreading  from  the  point  E  on  the 
edge  to  the  receiver  does  so  via  direct  and  ground-reflected  paths. 

The  ground  reflections  are  handled  by  simple  geometrical  acoustics; 
however,  the  "spherical  wave  reflection  coefficient"  is  used. 
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The  construction  of  the  "Edge-Plus-Images"  model  cen  best  be 
explained  with  the  help  of  Illustrations.  Figure  3.4  shows 
schematically  how  the  propagation  from  source  to  reciver  can  be 
modelled  as  a  superposition  of  four  distinct  half-plane  diffracted  ray 
paths.  In  Figure  3.4  (a),  the  "direct"  ray  from  source— to  edge— to 
receiver  does  not  interact  with  the  ground  and  so  is  represented  by 
the  standard  GTD  half-plane  equation.  The  diffraction  coefficient  can 
be  either  that  from  Equation  (3.6)  for  an  Impedance-covered  barrier  or 
that  from  Equation  (3.4)  for  a  rigid  barrier.  Is  the  equations  In 
Figure  3.4,  the  geometrical  spreading  factor  AO^.R^)  is  simply  (see 
Equation  [3.1]) 

A(Ra,Eb)  - 


In  Figure  3.4  (b),  the  source-ground-edge-receiver  ray  Is  still 
modelled  as  a  half-plane  diffraction  problem,  but  the  Incident  ray  is 
assumed  to  originate  from  the  source  "Image"  In  the  ground  Sj  and  Is 
correspondingly  modified  by  the  spherical  wave  reflection  coefficient 

^S^^1*Z1*R3^*  ln  th*  noct  *en*r*l  aanse,  this  reflection  coefficient 
Is  as  given  by  Equation  (2.140),  but  in  almost  all  practical 
situations,  the  first  term  ln  the  expansion  provides  sufficient 


accuracy.  Therefore,  the  model  has  been  Implemented  here  using  F^,  as 
defined  by  Equations  (2.149)  and  (2.150),  for  the  ground  reflections* 


As  Indicated,  the  factor  Qc 
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depends  on  the  reflection  angle  4^*  the 


ground  impedance  on  the  source  side  Z^,  and  the  image-to-edge  distance 


Figure  3.4  Decomposition  of  the  Overall  Problem  into  Four  Half-Plane  Diffraction 
Problems  in  Terms  of  the  Geometric  Theory  of  Diffraction. 

(b)  The  Reflected-Diffracted  Ray. 


Decomposition  of  the  Overall  Diffraction  Problem  into  Foi 
Problems  in  Terms  of  The  Geometric  Theory  of  Diffraction 
(c)  The  Dif fracted-Dlrect  Ray. 


Decomposition  of  the  Overall  Problem  into  Four  Half-Plane  Diffraction 
Problems  in  Terms  of  The  Geometric  Theory  of  Diffraction. 

(d)  The  Reflected-Diffracted-Reflected  Ray. 
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Rj.  It  is  interesting  to  note  that  the  edge-to-receiver  distance  R^ 
for  this  path  differs,  in  general,  from  the  analogous  distance  R^  for 
path  1.  This  is  a  consequence  of  the  law  of  edge  diffraction.  Since 
the  incident  ray  is  coming  from  the  source  image  and  not  the  source 
itself,  it  makes  a  different  angle  with  the  edge  and  defines  a 
different  "point  E"  along  the  edge.  Therefore,  the  exit  ray  to  the 
receiver  (making  an  equal  angle  with  the  edge)  must  correspondingly 
change  in  length.  The  diffraction  coefficients  reflect  this  fact 
also,  and  thus  the  notation  for  the  D^'s  has  subscripts  on  the  angle 
8* 

Similar  reasoning  applies  to  path  3,  the  source-edge-ground- 

receiver  ray,  where  QR  (ljj^Zj^g)  depends  on  the  properties  of  the 

ground  on  the  receiver  side  of  the  barrier.  To  model  the  final 

source-ground-edge-ground-receiver  path,  two  reflection  coefficients 

Qe  and  Q„  are  required,  as  indicated.  Note  that  since  a  different 
2  *2 

angle  6  is  again  defined  by  this  ray  path,  the  ground  reflection 
angles  and  are  generally  not  equal  to  ^  and  4^ »  respectively. 
The  total  diffracted  field  at  the  receiver  is 

4diff  "  *diff.  +  *diff,  +  $diff.  +  *diff .  .  (3.14) 

tot  12  3  4 

If  the  receiver  were  located  outside  of  the  acoustic  shadow  cast 
by  the  barrier-note:  each  ray  path  shown  in  Figure  3.4  defines  its 
own  shadow  boundaries— then  the  appropriate  non-dlffracted  geometrical 
fields  must  be  added  to  the  above  solution  to  obtain  the  total  field 
at  the  receiver.  That  is,  with  reference  to  Figure  3.3, 
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where  is  the  direct  field  from  source  to  receiver,  and  •  ^  ii 

the  barrier-reflected  field. 

Several  limiting  cases,  which  simplify  the  model  equations  in 
Figure  3.4,  may  be  of  practical  interest.  Thus, 


(a)  When  the  source  ground  is  rigid,  Q  -  Q  -1.0 

S1  S2 

(b)  When  the  receiver  ground  is  rigid,  Q  -  Q  -  1.0 

R1  R2 


(c)  When  the  source  is  on  the  ground,  Q 
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(d)  When  the  receiver  is  on  the  ground,  Q  -  Q  ; 

R1  R2 


r2  -  r6  -  r8  ; 


r5  and  R^  -  r7  ; 
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(e)  When  Che  source  and  receiver  are  boch  on  the  ground. 


$diff  -  0  $ 
a  2  diff. 
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diff 3  diff1 
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(f)  When  Che  source  and  receiver  are  boch  on  a  rigid  ground. 
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CHAPTER  IV 


NUMERICAL  RESULTS 

The  mathematical  solution  for  the  spherical  wave  propagation 
problem,  as  well  as  the  model  developed  for  barrier  diffraction,  are 
functions  of  several  geometrical  and  physical  parameters.  To 
determine  how  the  sound  field  at  the  receiver  predicted  by  the 
spherical  wave  theory  differs  from  that  predicted  from  the  classical 
plane  wave  theory,  numerical  results  for  several  representative 
source-receiver  geometries  will  be  presented  in  this  chapter.  For 
propagation  in  the  absence  of  a  barrier,  the  dependence  on  the 
parameter  kR,  the  reflection  angle  and  the  ground  Impedance  Z^, 
will  be  investigated.  When  a  diffracting  barrier  is  prese.'*  the 
dependence  of  the  field  on  barrier  height  and  di^  *,«  tion  angle  will 
also  be  studied. 

In  practice,  it  may  be  Important  to  know  when  a  problem  demands 
the  use  of  the  more  exact  spherical  wave  theory,  since  the  gain  in 
accuracy  could  be  significant.  Conversely,  the  plane  wave  solution  is 
extremely  simple,  and  its  accuracy  may  be  sufficient  enough  to  make 
its  use  economically  advisable.  This  is  particularly  true  for 
problems  in  which  a  continuous,  physically-extended  source  of  sound  is 
being  modelled  as  a  sequence  of  point  sources  (such  as  the  flow  of 
highway  traffic),  and  for  which  the  field  calculations,  therefore, 
must  be  performed  many  times. 
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4.1  Numerical  Results  for  the  Ground  Propagation  Problem 

This  section  presents  graphical  data  predicted  from  the  ground 
propagation  theories  given  in  Chapter  2.  The  ground  plane  is  assumed 
to  be  locally-reacting  with  its  Impedance  specified  by  real  and 
imaginary  parts  (since  the  time  convention  e~ia)t  has  been  adopted 
here,  a  positive  imaginary  part  X  of  the  impedance  ZN»R+iX  corresponds 
to  a  "springlike"  reactance).  Data  from  the  spherical  wave  theory 
derived  in  this  study  is  presented  in  two  forms,  representing  two 
levels  of  complexity.  First,  the  full  asymptotic  series  solution  is 
used  as  the  most  accurate  description  of  the  field;  second,  only  the 
first  term  in  this  series  is  employed,  as  an  approximation. 

Preliminary  tests  have  shown  that  both  of  these  forms  are  so  close  to 
the  exact  solution  for  many  practical  impedance  and  geometrical 
conditions  that  graphs  of  the  data  would  essentially  show  one  curve. 
For  this  reason,  the  values  of  the  parameters  chosen  for  many  of  the 
following  plots  are  deliberately  "atypical"  so  that  differences  in  the 
data  may  be  seen  and  studied. 

4.1.1  Dependence  on  the  Parameter  kR 

As  noted  in  the  derivation  in  Chapter  2,  the  asymptotic  series 
solution  for  the  point-to-point  ground  propagation  problem  approaches 
the  exact  (integral)  solution  as  the  parameter  kR^  increases  in  value. 
For  kR£  »  1  the  series — in  Equation  (2.130),  for  instance — is 
extremely  convergent  and  accurate  in  the  first  few  terms,  well  below 
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the  point  where  the  terms  begin  to  diverge.  As  kR^  gets  smaller,  more 
terms  are  necessary  for  accuracy,  yet,  correspondingly,  the  onset  of 
divergence  occurs  sooner.  Finally,  when  kR^  Is  very  small,  divergence 
Is  Immediate,  and  increasing  the  number  of  terms  only  degrades  the 
accuracy  of  the  solution. 

For  very  large  kR^ ,  the  curvature  of  the  wavefronts  becomes 
locally  ’’plane"  In  the  vicinity  of  the  ground  reflection  point,  and 
the  asymptotic  solution  should  approach  the  desired  plane  wave 
solution.  That  this  is  indeed  the  case  is  confirmed  by  the  numerical 
data. 

The  graphs  in  this  section,  as  well  as  most  of  the  others  in  this 

chapter,  compare  four  descriptions  of  the  sound  field  at  the  receiver. 

/ 

The  numerical  data  were  calculated  on  a  digital  computer  (IBM  Model 
3033  Processor)  from  (1)  a  fine-point  numerical  integration2*  of  the 
exact  Integral  expression  in  Equations  (2.93)  and  (2.101),  (li)  an 
"optimal"  number  of  terms  in  the  formal  asymptotic  series  solution 
given  by  Equation  (2.128),  (Hi)  the  first  term  in  the  latter  series 
as  represented  by  Equations  (2.149)  -  (2.151),  and  (iv)  the  plane-wave 
reflection  coefficient  form  of  the  solution.  Equations  (2.4)  and 
(2.5).  For  calculations  using  the  asymptotic  series,  a  truncation  of 
the  series  was  made  at  the  point  where  the  terms  just  begin  to 
diverge— although  always  including  a  minimum  of  two  or  a  maximum  of 
fifteen  terms.  The  solution  using  only  the  first  term  in  the  series 

24 

The  numerical  integration  was  performed  using  a  standard  five- 
point  Newton-Cotes  quadrature  technique  taking  a  finer  and  finer  mesh 
of  points  until  the  desired  degree  of  convergence  was  obtained. 


Is  referred  to  as  the  "F-term  solution**  in  the  following  discussion. 

To  facilitate  comparisons  of  the  data  within  each  graph,  the 
curves  in  this  chapter  have  been  drawn  as  straight  line  segments 
joining  discrete,  calculated  points.  Strictly  speaking,  then,  the 
values  on  a  line  segment  between  two  calculated  points  do  not 
represent  valid  data,  and  the  apparent  ** jaggedness **  of  the  curves  in 
some  regions  is  not  indicative  of  the  true  shape  of  the  data.  The 
number  of  calculated  data  points  may  vary  from  one  plot  to  the  next, 
but  they  are  always  equally-spaced  along  the  axis,  and  their  actual 
locations  can  usually  be  deduced,  if  necessary. 

The  data  for  the  case  where  both  the  source  height  (Sht)  and  the 
receiver  height  (Rht)  are  1.0  feet  above  the  impedance  plane,  and 
where  the  separation  (Sep)  is  38  feet  (corresponding  to  a  reflection 
angle  of  <|>-3°),  are  plotted  in  Figure  4.1  (a)  and  (b).  The  normalized 
ground  impedance  has  been  assigned  the  relatively  low  value  ■  (0.3, 
0.5).  The  vertical  axis  for  Figure  4.1  (a)  is  in  terms  of  the 
attenuation,  defined  as 

9 

Attenuation  ■  -  20  log^g  — ^t(^-  dB  ,  (4.1) 

air 

where  $tot  is  the  total  (direct  plus  reflected)  field  at  the  receiver, 

iVR 

and  $dir  ia  the  direct  field  only  (e  1/R^).  The  latter  is  the 
"free**  field  that  would  have  existed  at  the  receiver  in  the  absence  of 
the  ground.  The  same  quantity  defined  by  Equation  (4.1)  is  sometimes 
referred  to  as  the  "excess  attenuation  relative  to  spherical 


Figure  4.1  Calculated  Data  for  (a)  Attenuation,  and 

(b)  Reflected/Direct  Ratio  as  a  Function  of  kR,.  Exact, 
Asymptotic  Series,  F-term,  and  Plane  Wave  Solutions. 
Small  Reflection  Angle,  Low  Impedance. 


spreading. "  Finally,  it  should  be  noted  that  a  negative  value  of 
attenuation  means  that  the  received  field  has  been  enhanced  relative 
to  free  field  level. 
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As  expected,  and  as  is  readily  apparent  from  either  Figure  4.1 
(a)  or  (b),  the  plane  wave  solution  (long-dashed  line)  is  not  very 
accurate  for  this  near-grazing  incidence  geometry.  On  the  other  hand, 


the  asymptotic  series  solution  (solid  line)  agrees  well  with  the  exact 
integral  (short-dashed  line)  for  IcR^  >3,  a  surprisingly  low  value  for 
kRj.  Below  this  value,  the  series  is  divergent,  and  the  predictions 
become  erratic.  The  medium-dashed  line  shows  that  using  only  the 
first  term  in  the  full  solution— the  F-terrn  solution—  gives  very  good 
agreement,  differing  by  a  few  dB  for  medium  values  of  kR^.  It  is 
remarkable,  though,  that  the  F-term  solution  yields  very  good 
predictions  down  to  kRg  ■  0.1,  corresponding,  for  the  distances 
considered  here,  to  a  frequency  of  0.03  Hz. 

All  four  solutions  eventually  agree  when  IcR^  is  sufficiently 
large.  Physical  Intuition  predicts  this— the  wavefronts  "look"  like 
plane  waves— but  it  can  also  be  deduced  from  the  governing  equations. 
For  example,  with  reference  to  Equation  (4.4),  as  kR^  gets  large,  \ 
also  gets  large,  and  F^  can  be  shown  to  have  the  asymptotic  expansion 
(see,  in  part,  Abramowltz  &  Stegun,  1964,  section  7.1.23): 

.2  1  1*3  1*3*5 

F  *  2±&\  e  H(-Io(X)  ] - -  (4.6) 

2X2  U4  8*6 

where  the  Heaviside  operator 

* 

1  ;  Im(z)  >  0 

H(  z  ]  ■  1  (4.7) 

0  ;  Otherwise 

is  necessary  to  account  for  the  fact  that  erfc(-z)  »  2  -  erfc(z). 

Thus,  Fj  becomes  negligible  as  kR2  gets  large,  and  Equation  (4.4) 
reverts  to  the  plane  wave  solution. 
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The  attenuation  plotted  In  Figure  4.2  (a)  and  the  ratio  of  the 
fields  In  Figure  4.2  (b)  are  for  the  sane  source-receiver  geometry  as 
In  the  previous  case,  but  for  a  higher  ground  Impedance,  Z„  -  (2.0, 
3.0).  Two  observations  can  immediately  be  made  when  Figure  4.1  and 
Figure  4.2  are  compared.  First,  the  exact,  the  full  asymptotic 
series,  and  the  first-term-only  solutions  are  in  very  close  agreement, 
even  for  values  of  kR2  down  to  0.1.  In  fact,  for  k&2  >,4,  the 
predictions  are  essentially  identical.  Second,  the  "overall*'  level  of 
attenuation  is  lower  for  the  higher  Impedance  ground.  This  trend  is 
also  not  surprising,  and  in  the  limit  of  a  rigid  ground,  both  the 
classical  theory  and  the  present  spherical  wave  theory  predict  a 
"gain"  of  -6  dB. 

Some  further,  less  obvious,  conclusions  can  be  drawn  from  the 
results  of  Figures  4.1  and  4.2.  The  peak  of  maximum  attenuation 
(centered  at  about  kR2«20  for  ZN»[0.3,0.5J  and  at  about  kR2»150  for 
Zjj-[2. 0,3.0])  is  not  a  result  of  propagation  path  length  Interference, 
but  rather  Is  strictly  an  "impedance  effect."  For  a  fixed  geometry, 
this  maximum  will  move  to  higher  and  higher  frequencies  as  the 
Impedance  increases  and  will  eventually  disappear  or  "merge"  with  the 
first  true  path  length  interference  peak  In  the  limit  when  the  ground 
becomes  rigid. 

Figure  4.3  illustrates  the  numerical  results  for  the  case  where 
source  and  receiver  are  both  directly  on  the  ground  (  ^  »  0°)  and 
separated  by  57  feet.  The  impedance  is  the  same  as  in  Figure  4.1, 

■  (0.3,  0.5).  The  plane  wave  curve  is  not  shown  here,  since  the 


predicted  attenuation  would  be  infinite  for  all  values  of  kR^ 

(Rp  ■  -1,  Rj-R^.  The  asymptotic  solution  again  gives  extreme 
accuracy  for  kR^  ^  4,  while  the  F-term  solution  differs  by  only  a  few 
dB  over  most  of  the  range. 

An  interesting  fact  can  be  observed  in  either  Figure  4.1  (a)  or 
in  Figure  4.3.  That  is,  for  very  low  kR2,  the  exact  solution  (as  well 
as  the  F-term  solution)  predicts  a  gain  of  the  order  -7  dB.  Thus,  the 
field  at  the  receiver  is  actually  higher  than  that  which  would  exist 
if  the  ground  were  perfectly  reflective.  This  is  again  a  direct 
effect  of  the  finite  ground  impedance  and  has  been  referred  to  in 
various  ways  in  the  literature,  among  them  the  “surface  wave  effect,** 
a  “ducting  of  energy, “  a  “focallzatlon”  phenomenon,  and  a  “surface 
layer  effect."  Whatever  the  terminology,  this  phenomenon — the  gain 
could  be  higher  than  7  db—is  a  testament  to  the  fact  that  the 
reflected  wave  is  not  spherically  symmetric,  but  that  Instead  energy 
is  "re-radiated"  along  the  entire  boundary.  This  re-radiation  is  with 
different  phases  from  different  points  along  the  surface  and  could 
result  in  constructive  interference  at  certain  field  points. 

Some  conclusions  about  the  effect  of  Increasing  the  reflection 
angle  can  be  drawn  from  the  next  pair  of  plots  in  Figure  4.4  (a)  and 
(b).  Here,  the  source  and  receiver  have  been  raised  to  five  feet 
above  the  ground  OP-IO0),  and  the  impedance  is  ZN-(0.3,  0.5)  in  the 
former  and  Z^»(2.0,  3.0)  in  the  latter.  The  counterpart  cases  are 
those  in  Figures  4.1  (a)  and  4.2  (a).  Comparing  Figure  4.1  (a)  to 
Figure  4.4  (a)  and  Figure  4.2  (a)  to  Figure  4.4  (b)  shows  that  the 
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peak  of  maximum  attenutlon  shifts  to  lower  values  of  kR^  as  ’P 
Increases.  Inherent  in  this  trend  is  the  fact  that  the  range  of  kR^ 
for  which  the  plane  wave  solution  becomes  valid  begins  at  a  lower 
value  when  the  source  and  receiver  are  situated  at  a  higher  elevation 
above  the  ground.  Furthermore,  for  each  impedance  case,  the  overall 
attenuation  decreases  as  gets  larger,  due  to  the  fact  that  the 
ground  appears  ’’harder"  at  larger  reflection  angles.  It  might  be 
Interesting  to  note  that  the  second  peak  in  both  Figure  4.4  (a)  and 
(b)  (at  kR^  3  500)  is,  in  fact,  due  to  path  length  differences  and 
hence  does  not  shift  significantly  as  \p  goes  from  3°  to  10°  (K, 

=  R,j) .  Conversely,  the  impedance  peak  shifts  by  more  than  an  order 
of  magnitude. 

To  complete  this  section,  numerical  results  for  a  case  of  a  very 
large  angle  of  reflection  and  a  very  low  impedance  are  presented  in 
Figure  4.5  (a)  and  (b).  The  source  height  is  10  feet;  the  receiver 
height  is  20  feet;  the  separation  is  30  feet  (^45°);  and  the  ground 
Impedance  is  ZN-(0.4,  0.3).  Here,  the  lmpedance~effect  attenuation 
peak  has  shifted  to  a  very  low  value  of  kR2,  and  path  length 
Interference  prevails  for  kR^  >  10.  The  plane  wave  solution  yields 
predictions  as  accurate  as  the  F-term  predictions  for  all  kR^  greater 
than  about  5. 

Table  4.1  summarizes  the  geometry  and  impedance  conditions  for 
the  cases  plotted  in  this  section. 
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Figure  4.5  Calcualted  Data  for  (a)  Attenuation  and  (b)  Reflected/Direct 
Ratio  as  a  Function  of  kR2-  Exact,  Asymptotic  Series,  F-term,  and  Plane 
Wave  Solutions.  Large  Reflection  Angle,  Low  Impedance. 


Table  4.1 


4.1.2  Dependence  on  Reflection  Angle 


To  investigate  the  dependence  of  reflection  angle  ip  on  the 
received  acoustic  field,  it  is  necessary  to  hold  the  parameter  klL, 
constant.  Thus,  the  approach  illustrated  in  Figure  4.6  was  used  for 
generating  the  numerical  data  plotted  in  this  section.  The  source 
image-to-receiver  distance  R2,  the  source  height,  and  the  frequency 
are  held  constant  as  the  receiver  point  revolves  around  the  image.  In 
this  way,  ip  varies  from  its  initial  position  with  the  receiver  on  the 
ground  to  its  final  value  of  90°.  Of  course,  R^  also  will  vary  with 
the  angle  ,  but  all  of  the  parameters  in  the  reflected  field  terms 
will  remain  constant. 

The  results  plotted  in  Figure  4.7  and  Figure  4.8  compare  two 
values  of  the  parameter  kR^  for  a  low  value  of  ground  impedance  Here, 
the  source  height  is  fixed  slightly  above  the  ground  at  0.1  feet;  R2 
is  12  feet;  -  (1.2,  1.8);  and  kR2  ■  6.7  (frequency  ■  100  Hz.)  in 
Figure  4.7  and  kR2  ■  33.5  (frequency  ■  500  Hz.)  in  Figure  4.8. 

The  asymptotic  series  solution  is  indistinguishable  from  the 
exact  integral  solution  for  all  values  of  ip,  and  only  very  slight 
differences  can  be  detected  with  the  F-term  form  of  solution.  On  the 
contrary,  the  plane  wave  solution  is  totally  inaccurate  for  small 
grazing  angles  yet  improves  as  ip  increases.  For  the  low  value  of  kR2 
■  6.7,  this  improvement  is  not  significant  until  ip  is  greater  than 
30°,  whereas  for  the  modestly-high  value  of  kR2  -  33.5,  the  plane  wave 
solution  is  highly  accurate  for  values  above  i|>*5° 
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kR2  *  6.7 

Zn*U.2.  1.8) 

-  Exact 

-  Asym.  Series 

- F-term 

- Plane  Wave 
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Reflection  Angi®,  D«gr««s 

Figure  4.7  Calculated  Data  for  (a)  Attenuation  and  (b)  Reflected/Direct 
Ratio  as  a  Function  of  Reflection  Angle.  Exact,  Asymptotic  Series 
F-term,  and  Plane  Wave  Solutions.  Small  Idt^,  Low  Impedance. 
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The  response  curves  of  attenuation  vs.  reflection  angle  will  be 
different  for  different  ground  impedances.  In  Figure  4.9  are  plotted 
four  sets  of  curves  corresponding  to  four  values  of  impedance  that 
differ  by  factors  of  ten  from  each  other.  The  plane  wave  solution  is 
compared  to  the  F-term  solution  only,  since  the  latter  has  been 
consistently  close  to  the  exact  and  asymptotic  solutions.  The  value 
of  kR^  has  been  fixed  deliberately  low  to  force  differences  in  the 
solutions  so  that  the  trends  in  the  data  could  be  observed.  Here,  the 
source  height  is  at  1.0  feet;  is  fixed  at  12  feet;  kR^  ■  3.35 
(frequency  ■  50  Hz.);  and  the  four  values  of  impedance  are  as  shown  on 
the  graph.  As  expected,  the  bottom  two  sets  of  curves  show  that  as 
the  Impedance  increases  from  Z^B(1.0,  0.5  )  to  ZN-(10.0,  5.0),  the  two 
solutions  fall  closer  into  accord.  However,  the  top  two  sets  of 
curves  show  the  opposite  trend.  For  the  extremely  low  impedance 
ZN-(0. 01 ,0.005) ,  the  solutions  are  identical,  yet  as  the  Impedance 
Increases  tenfold  to  ZN-(0.1,0.05),  the  predictions  begin  to  differ. 
This  is  an  interesting  result,  but  it  should  not  be  surprising.  The 
plane  wave  theory,  based  on  the  image  method,  is  exact  for  either 
infinitely  rigid  or  pressure  release  surfaces.  Therefore,  as  Z-+-0, 
the  plane  wave  solution  becomes  more  accurate,  and  it  should  agree 
with  the  spherical  wave  F-term  solution. 

The  geometry  and  impedance  conditions  for  the  cases  plotted  in 


this  section  are  summarized  in  Table  4.2. 


Figure  4.9  Comparison  of  Calculated  Data  for  Four  Values  of  Ground  Impedance,  as  a  Function  of 
Reflection  Angle.  F-term  and  Plane  Wave  Solutions. 


4.1.3  Received  Spectra  Over  a  Practical  Ground  Cover 


In  several  of  the  previous  graphs,  the  value  of  ground  Impedance 
was  held  constant  as  kR^  or  frequency  varied.  While  this  may  be 
necessary  for  studying  the  sensitivity  of  various  parameters,  the 
impedance  of  actual  surfaces  is  rarely  frequency-independent.  The 
problem  in  practice  is  more  often:  given  a  ground  surface  between  a 
source  and  a  receiver  characterized  by  its  impedance  measured  at 
several  discrete  frequencies,  predict  the  attenuation  at  the  receiver. 
A  common  ground  cover  is  outdoor  grass,  and  one  particular  type — the 
so-called  NRC  grass— has  the  third-octave  values  of  impedance  shown 
in  Table  4.3.  Using  the  values  of  frequency  and  impedance  given  in 
the  table,  plots  of  attenuation  versus  frequency  have  been  generated 
and  are  plotted  in  Figures  4.10  -  4.12.  The  first  two  of  these  graphs 
shows  the  effect  of  symmetrically  raising  both  source  and  receiver, 
while  keeping  the  horizontal  separation  constant.  Thus,  with  a 
separation  of  50  feet,  the  source  and  receiver  move  from  0.1  feet 
above  the  ground  (^«0.23°)  in  Figure  4.10,  to  1.0  feet  (ty-2.3°)  in 
Figure  4.11  (a),  and  to  2.0  feet  (ty»4.5°)  in  Figure  4.11  (b). 

Although  the  impedance  is  not  constant  here,  the  trend  of  the 
attenuation  maximum  moving  to  lower  frequencies  and  the  overall 
attenuation  decreasing  with  increasing  <l>  noted  in  section  4.1.2  is 
still  apparent.  Again,  the  plane  wave  solution  becomes  more  accurate 

^  The  impedance  of  this  outdoor  surface  was  measured  by  several 
methods  at  the  National  Research  Council  in  Ottawa  (Embleton,  et  al., 
1975);  hence  the  abbreviation  in  the  name. 
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Figure  4.10  Calculated  Data  for  (a)  Attenuation  and  (b)  Reflected/ 
Direct  Ratio  as  a  Function  of  Frequency  for  Propagation 
Over  Grass.  F-term  and  Plane  Wave  Solutions. 

Sht-0.1,  Rht-0.1,  Sep-50.0,  Reflection  Angle-0. 23°. 
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Figure  4.12  Calculated  Data  for  (a)  Attenuation  and  (b)  Reflected/ 
Direct  Ratio  as  a  Function  of  Frequency  for  Propagation 
over  Grass.  F-tero  and  Plane  Wave  Solutions. 

Sht-1.0,  Rht-1.0,  Sep-500.0,  Reflection  Angle-  0.23°. 


for  larger  angles  of  reflection.  Of  practical  Interest  Is  the  fact 
that  the  grass-covered  ground  appears  rigid  (attenuation  -  -6  dB)  for 
frequencies  below  about  800  Hz.  The  low-frequency  noise  of 
automobiles  and  trucks  is  not  well-attenuated  over  grassland. 

The  data  in  Figure  4.12  are  for  source  and  receiver  both  1.0  feet 
above  the  ground  and  separated  by  500  feet.  In  this  case,  then,  the 
reflection  angle  is  the  same  as  for  the  case  plotted  in  Figure  4.10, 
namely  if>  -0.23°,  but  the  distance  is  larger.  With  reference  to 

Figure  4.10  (b)  and  4.12  (b),  it  can  be  seen  that  the  reflected  field 
for  the  plane  wave  has  not  changed  at  all  since  Rp  is  independent  of 
R2;  however,  the  spherical-wave  reflected  field  is  quite  different  in 
Figure  4.12  than  in  Figure  4.10  because  the  F-term  itself  depends 
explicitly  on  R£  (see  Equation  4.4). 

To  conclude  this  section,  data  for  near-grazing  propagation  over 
a  ground  surface  comprising  an  indoor-outdoor  carpet  material  is 
presented  in  Figure  4.13.  The  impedance  values  for  the  material  are 
given  in  Table  4.4.  This  particular  carpet  has  found  use  in  reduced- 
scale  model  experiments  (Lawther,  et  al.,  1980),  and  hence  the 
frequencies  in  both  Table  4.4  and  in  Figure  4.13  are  higher  than  those 
presented  thus  far.  The  interesting  feature  in  the  plotted  data  is 
that  the  "gain”  is  well  in  excess  of  -6  dB  for  several  values  of 
frequency.  This  is  most  clearly  shown  in  Figure  4.13  (b)  where  the 
magnitude  of  the  reflected  field,  given  by  Equation  4.4,  can  be  more 
than  3  times  as  great  as  the  direct  field. 
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Figure  4.13  Calculated  Data  for  (a)  Attenuation  and  (b)  Reflected/ 
Direct  Ratio  as  a  Function  of  Frequency  for  Propagation  over  Carpet. 
F-tena  and  Plane  Wave  Solutions.  Sht»0.33,  Rht-0.25,  Sep-24.0,  <|>-1.4°. 


Table  4.4 


Real  and  Imaginary  Parts  of  the  Impedance 
of  Indoor-Outdoor  Carpet  Material 
for  Third-Octave  Values  of  Frequency. 


Frequency 


Imaginar 
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The  geometry  and  Impedance  conditions  for  the  plots  given  in  this 
section  are  summarized  in  Table  4.5 

4.1.4  Results  for  Perpendicular  Incidence 

The  exact  solution  to  the  problem  of  point  source  reflection  by 
an  impedance  plane  for  perpendicular  incidence  is  given  in  Section  2.4 
in  terms  of  the  exponential  integral.  For  small  arguments  the 
exponential  Integral  could  be  expressed  in  a  convergent  series 
expansion,  and  for  large  arguments,  it  can  be  written  in  terms  of  its 
asymptotic  expansion.  As  has  been  shown,  the  first  term  in  the 
asymptotic  series  will  generate  the  plane  wave  solution  given  by 
Equation  (2.159). 

Plotted  in  Figure  4.14  are  the  magnitudes  of  the  ratios,  in  dB, 
of  (a)  the  total  field  at  the  receiver  to  the  free  field  and  (b)  the 
reflected  field  at  the  receiver  to  the  free  field  for  a  source  3.0 
feet,  and  a  receiver  1.0  feet,  above  the  ground  plane,  respectively. 
The  impedance  of  the  ground  is  ZN  -  (1.4,  1.3).  The  short-dashed 
curve  represents  the  exact  solution  to  the  perpendicular  incidence 
problem,  obtained  by  numerically  integrating  Equation  (2.152).  The 
medium-dashed  line  is  the  calculated  solution  using  the  first  two 
terms  in  the  convergent  series  expansion  of  the  exponential  integral, 
Equation  (2.156),  and  the  long-dashed  line  is  the  solution  using  the 
first  two  terms  in  the  asymptotic  series,  Equation  (2.157).  Finally, 
the  solid  line  curve  is  the  plane  wave  solution. 
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It  can  be  seen  that  the  predictions  using  the  first  two  terms  in 
the  convergent  series  is  very  accurate  for  kR^  <  1  but  veers 
erratically  for  values  above  that.  When  the  first  two  terms  in  the 
asymptotic  series  are  used,  on  the  other  hand,  excellent  agreement 
with  the  exact  solution  is  obtained  for  kR^  ^3,  whereas  the  plane 
wave  solution  (referring  to  Figure  4.14  [b])  requires  kR2  ^50,  or  so. 
It  should  be  emphasized  that  these  results  are  for  a  relatively  low 
value  of  Impedance,  and  that  for  most  practical  impedance  surfaces  (in 
acoustical  studies  in  air)  the  plane  wave  solution  will  be  much  more 
in  accord  with  the  exact  solution. 

Whether  the  value  of  kRj  is  small  or  large,  there  is  one  ground 
impedance  region  for  which  the  plane  wave  solution,  if  used,  should  be 

applied  with  caution.  That  is  where  the  value  of  admittance  8 - ►  1. 

From  Equation  (2.159),  it  can  be  seen  that  the  reflected  field 
vanishes  for  B~l,  and  the  attenuation  becomes  zero.  Physically,  this 
is  simply  the  ideal  "matched"  condition  where  Z^-l  (unnormalized 
Impedance  -  pc),  and  no  energy  is  reflected.  The  spherical  wave 
theory,  however,  does  predict  some  reflected  energy  for  this  case. 
Figure  4.15  (a)  and  (b).  show  the  familiar  "standing  wave  pattern" 
that  is  obtained  when  the  source  is  held  fixed  at  some  distant  point 
above  the  plane  (here,  100  feet)  and  the  receiver  is  moved  away  from 
the  surface.  Experiments  based  on  this  standing  wave  pattern  (using  a 
pressure  microphone  or  hydrophone)  are  often  conducted  to  calculate 
the  value  of  surface  Impedance.  The  ground  impedance  in  Figure  4.15 
(a)  is  ZN  -  (1.0,  0.05)  and  in  Figure  4.15  (b)  is  ZM  -  (1.0,  0.005). 


Although  the  scales  on  the  vertical  axes  are  greatly  expanded,  the 
difference  In.  the  theoretical  predictions  could  be  important  in 
practice. 

4.2  Results  for  the  Barrier  Diffraction  Problem 

This  section  presents  numerical  results  from  the  combined  sound 
propagation-barrier  diffraction  model  described  in  Chapter  3.  The 
ground  on  either  side  of  the  barrier  is  again  assumed  to  be  locally¬ 
reacting,  and  an  arbitrary  impedance  may  be  assigned  to  each  side 
separately.  The  source  and  receiver  can  be  located  at  any  non- 
coincident  points,  although  the  GTD  method  may  become  inaccurate  if 

either  is  closer  than  a  few  wavelengths  from  the  edge.  The  barrier 

26 

itself  will  be  assumed  rigid  for  the  present  calculations.  The 
relevant  parameters  to  be  investigated  in  this  section  are  again  kR 
and  ground  impedance  Z^,  but  in  addition,  the  dependence  of  the 
predictions  on  the  barrier  height  and  the  diffraction  angle  $  will  be 
examined . 

The  attenuation  as  defined  in  the  previous  section  will  again  be 
plotted,  since  it  is  a  useful  quantity  for  comparing  one  solution  to 
another  or  for  comparing  one  geometry  to  another.  In  noise  control 

26 

Although  the  theoretical  solution  for  an  impedance-covered 
half-plane  was  presented  in  Chapter  3,  prior  sensitivity  studies 
(Hayek,  Lawther,  Kendig,  &  Slmowitz,  1978)  have  shown  that  the 
differences  in  the  sound  field  for  the  impedance  barrier  relative  to 
the  rigid  barrier  were  not  significant  for  most  practical  cases. 
Computations  using  the  diffraction  coefficient  given  by  Equation 
(3.6),  then,  were  deemed  unecessary  for  the  purposes  of  the  present 


practice,  however,  a  more  informative  quantity  is  the  "insertion  loss’* 
of  the  barrier,  which  relates  the  actual  field  at  the  receiver  to  the 
field  that  would  have  existed  there  in  the  absence  of  the  barrier  (but 
in  the  presence  of  the  ground).  Thus,  graphs  of  the  data  will  also  be 
given  here  in  terms  of  the  insertion  loss,  defined  formally  as 

barrier 

Insertion  Loss  ■  -20  log^Q  ^ -  dB  ,  (4.8) 

groundtot 


where  $karrier  is  the  total  field  at  the  receiver  point  when  the 
rtot 


barrier  is  present  and  $ _ .  is  the  received  field  over  the 

ground tQt 


impedance  plane  alone.  Note  that  the  latter  field  must  be  calculated 
separately,  using  the  ground  propagation  theory  of  Chapter  2.  Because 
the  Edge-Plus-Images  model  uses  only  the  first  term  in  the  asymptotic 
series  solution  to  compute  the  ground  reflections,  the  calculation  of 


$ around  for  the  insertion  loss  data  has  also  been  performed  with 

8  tot 


only  the  first  term  (the  solution  given  in  Section  2. 3. 5 .5). 

Some  interesting  features  of  the  Edge-Plus-Images  model  should 
perhaps  be  pointed  out  before  discussing  the  numerical  results. 

First,  because  the  critical  angles  defining  the  acoustic  diffraction 
are  measured  from  the  top  of  the  barrier  edge  (refer  to  Figure  3.2), 
the  actual  (absolute)  location  of  the  ground  is  of  secondary 
Importance.  That  is,  from  the  standpoint  of  the  model,  the  ground 
need  not  be  horizontal.  The  configuration  depicted  in  Figure 
4.16  (a),  which  might  represent  a  barrier  on  a  hillside,  is 
mathematically  identical  to  that  shown  in  Figure  4.16  (b),  which  could 
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represent  a  “tilted”  barrier  on  a  horizontal  ground.  Furthermore, 

since  the  ground  may  have  different  “slopes”  on  either  side  of  the 

barrier,  geometries  such  as  that  in  Figure  4.16  (c)  can  be  addressed 

by  the  model  (provided  the  source  and  receiver  “images”  in  the  ground 

remain  on  the  same  side  of  the  extended  half-plane  as  the  source  and 

receiver  themselves).  Finally,  since  the  location  of  the  edge  point  E 

(Figure  3.1)  determines  the  diffraction  angles,  the  “height”  of  the 

barrier  need  not  be  constant  relative  to  the  ground.  Therefore,  a 

very  general  configuration  such  as  that  illustrated  in  Figure  4.16  (d) 

27 

can  be  solved  by  the  Edge-Plus-Images  model. 


4.2.1  Dependence  on  the  Parameter  kR 


Figure  4.17  (a),  (b),  and  (c)  shows  the  geometry  of  the  barrier 
problem  and  defines  the  parameters  which  will  appear  on  many  of  the 
graphs  in  this  section.  Rather  than  investigate  the  dependence  on  kR^ 
and  kR^  separately,  the  following  data  will  be  plotted  in  terms  of  kR, 
where  R  is  defined  as  R  ■  R^  +  R^.  Furthermore,  in  all  of  the 
following  cases,  the  distances  R^  and  have  been  kept  approximately 
equal  in  value.  The  data  plotted  in  Figure  4.18  (a)  and  (b)  is  for  a 
source  located  0.2  feet  above  the  ground  (Sht)  and  "offset"  from  the 
barrier  by  12.0  feet  (Soff),  while  the  receiver  is  located  0.1  feet 


27 

Of  course,  some  basic  solid  and  analytic  geometry  is  necessary 
for  computing  the  locations  in  space  of  the  various  image  points  and 
propagation  paths.  Included  in  the  computer  program  implementation  of 
the  present  model  are  several  geometrical  subroutines  for  performing 
these  tasks. 
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Figure  4.18  Calculated  Data  for  (a)  Attenuation  and  (b)  Insertion  Loss 
as  a  Function  of  kR.  Q-term,  Plane  Wave,  and  Half-Plane  Solu¬ 
tions.  Small  Source  and  Receiver  Ht.,  Small  Bht.,  Low  Z„. 


above  the  ground  (Rht)  and  offset  by  12.0  feet  (Roff).  The  barrier 
height  (Bht)  is  1.0  feet  and  the  ground  impedance  is  ZN  =  (1.5,  2.0) 
on  each  side.  As  shall  be  true  for  most  of  the  plots  in  this  section, 
the  source-receiver  ray  is  perpendicular  to  the  barrier  (zQ  =  0  and 
8  -  90°). 

Three  curves  are  drawn  corresponding  to  calculated  data  using  (i) 
the  "half-plane"  theory  (long-dashed  line),  (ii)  the  "plane  wave” 
theory  (short-dashed  line),  and  the  (iii)  "Q-term"  theory  (solid 
line).  For  the  half-plane  computations,  only  the  "direct”  diffraction 
path  (Path  1  in  Figure  3.4  [a])  is  used,  as  though  the  ground  were 
absent  and  a  semi-infinite  barrier  existed.  The  plane  wave  solution 
accounts  for  all  the  diffraction  paths,  but  the  plane  wave  reflection 
coefficient  given  by  Equation  (2.4)  is  used  in  place  of  the  spherical 
wave  reflection  cofflclent  Q.  The  data  computed  with  the  Q-term 
solution  is  the  most  accurate  of  the  three  and  employs  the  first  term 
in  the  spherical  wave  reflection  coefficient  given  by  Equation  (2.40) 
(that  is,  where  "F"  is  simply  replaced  by  "F^  from  Equation  [2.150]). 
In  all  cases,  the  diffraction  coefficient  Dg  is  that  given  by  Equation 
(3.4). 

The  curves  show  some  interesting  features.  The  Q-term  shows  a 
broad  peak  in  the  attenuation  (and  insertion  loss)  similar  to  that 
found  in  the  previous  ground  propagation  graphs.  This  is  again  due  to 
an  "impedance  effect,"  as  path  length  interference  can  not  occur  in 
this  kR  region.  For  very  large  kR,  the  plane  wave  solution  agrees 
with  the  Q-term  solution,  as  expected,  since  in  this  region  the 
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spherical  wave  reflection  coefficient  reduces  to  the  plane  wave 
reflection  coefficient.  However,  the  two  predictions  are  very 
different  for  smaller  values  of  kR.  This  can  be  explained  as  follows. 
Due  to  the  short  barrier  height,  the  ground  reflection  angles  on 
either  the  source  or  receiver  sides  are  small  (approximately  5°),  and, 
consequently,  the  plane  wave  reflection  coefficients  approach  -1.0. 
Therefore,  instead  of  showing  a  gain  over  free  field  as  predicted  by 
the  Q-term  solution,  use  of  the  plane  wave  reflection  coefficient 
shows  an  attenuation  of  about  20  d&. 

The  curve  generated  from  the  half-plane  theory  is  inaccurate  over 
the  full  range  of  kR  considered  here.  The  only  discernible  trend  is 
the  gradual  increase  in  attenuation,  reflecting  the  fact  that  the 
diffraction  coefficient  in  Equation  3.4  decreases  with  increasing  kR. 

The  insertion  loss  data  in  Figure  4.18  (b)  reveals  an  interesting 

fact.  Namely,  for  large  kR  the  predicted  insertion  loss  is  negative, 

indicating  that  the  field  at  the  receiver  behind  the  barrier  is 

28 

actually  greater  that  it  would  be  without  the  barrier.  The 
explanation  for  this  phenomenon  is  simple:  in  the  absence  of  the 
barrier,  the  reflection  angle  from  source  to  receiver  is  of  the  order 
of  0.7°,  and  the  "propagation  loss"  over  the  relatively  "soft" 
impedance  plane  would  be  large.  The  barrier  essentially  forces  the 
reflection  angles  higher  (about  5°),  causing  the  ground  to  appear 

28 

This  is  also  true  for  low  values  of  kR,  but  the  attenuation  is 
likewise  predicted  to  be  negative  there.  Moreover,  both  the  ground 
propagation  theory  and  the  GTD  are  extremely  accurate  for  "large  kR," 
and  hence  the  assertion  of  negative  insertion  loss  can  be  made 
confidently. 
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"harder,”  and  thereby  reducing  the  propagation  loss. 

The  numerical  data  In  Figure  4.19  (a)  and  (b)  Is  for  the  same 
geometry  as  In  the  previous  case,  only  here  the  ground  Impedance  has 
Increased  by  a  factor  of  10.  Immediately  apparent  is  the  fact  that 
the  overall  attenuation  is  lower  for  the  harder  ground.  Actually,  as 
was  the  trend  in  the  ground  propagation  studies  in  the  previous 
section,  the  attenuation  peak  has  shifted  to  higher  values  of  kR. 

Also  apparent  is  that  the  plane  wave  predictions  are  closer  to  the  Q- 
term  predictions  (maximum  deviation  5  dB,  as  opposed  to  28  dB  in 
Figure  4.18).  This  is  not  surprising  since  the  influence  of  ground 
absorption  diminishes  as  the  ground  impedance  increases.  Again,  if 
the  ground  impedance  were  increased  continuously,  the  two  solutions 
would  eventually  give  identical  results. 

Unlike  the  case  for  the  soft  ground,  the  insertion  loss  in  Figure 
4.19  (b)  is  appreciable  at  high  frequencies.  This  is  because  the  loss 
due  to  simple  propagation  would  not  be  significant  over  this  ground, 
so  the  "shadowing”  effects  of  the  barrier  appear  stronger  by 
comparison. 

The  attenuation  and  insertion  loss  curves  presented  next  in 
Figures  4.20  and  4.21  show  some  preliminary  effects  of  increasing  the 
barrier  height  (this  dependence  will  be  discussed  in  detail  in  the 
next  section).  All  parameters  are  the  same  as  those  in  Figures  4.18 
and  4.19,  respectively,  only  the  barrier  height  is  now  8.0  feet 
instead  of  1.0  feet.  For  either  impedance  considered,  the  range  of  kR 
for  which  the  plane  wave  solution  is  in  accord  with  the  Q-term 


Sht  =  0.2  ft 
Rht  =  0.1  ft 
Soft  =  12.0  ft 
Roff  =  12.0  ft 
Bht  =  1.0  ft 
ZN  =  (15.0,  20.0) 


Q-term 
Plane  Wave 
Half-Plane 


Sht  =  0.2  ft 

Rht  =  0.1  ft 

Soft  =  12.0  ft 

Roff  -  12.0  ft 

Bht  =  1.0  ft 

ZN  *  (15.0,  20.01 

/  / 

-  Q-term 

- Plane  Wave 

- Half-Ptane 

Figure  4.19  Calculated  Data  for  (a)  Attenuation  and  (b)  Insetion 

Loss  as  a  Function  of  kR.  Q-term,  Plane  Wave,  and  Half 
Plane  Solutions.  Small  Source  and  Receiver  Heights, 
Small  Barrier  Height,  High  Impedance. 


Sht  =  0.2  ft 
Rht  =  0.1  ft 
Soft  =  12.0  ft 
Roff  =  12.0  ft 
Sht  -  3.0  ft 
Z,,  *  (1.5,  2.0) 


— —  Q-term 

- Plane  Wave 

- Half-Plane 


\ 

\ 

\ 


—  Q-term 

—  Plane  Wave 

—  Half-Plane 


Sht  -  0.2  ft 
Rht  =  0.1  ft 
Soff  =  12.0  ft 
Roff  =  12.0  ft 
Sht  =  8.0  ft 
ZN  =  (1.5.  2.0) 


\  V 
\  / 


Figure  4.20  Calculated  Data  for  (a)  Attenuation  and  (b)  Insertion 

Loss  as  a  Function  of  kR.  Q-term,  Plane  Wave,  and  Half- 
Plane  Solutions.  Small  Source  and  Receiver  Heights, 
Large  Barrier  Height,  Low  Impedance. 
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Figure  4.21  Calculated  Data  for  (a)  Attenutation  and  (b)  Insertion 

Loss  as  a  Function  of  kR.  Q-term,  Plane  Wave,  and  Half- 
Plane  Solutions.  Small  Source  and  Receiver  Heights, 
Large  Barrier  Height,  High  Impedance. 


solution  begins  at  a  lower  value  for  the  higher  barrier.  That  is,  for 
ZN  -  (1.5,  2.0),  the  1-foot  barrier  dictates  kR  _>  700  for  agreement, 
whereas  the  8-foot  barrier  requires  kR  >,  50.  For  -  (15.0,  20.0),  a 
large  kR  is  still  necessary  for  the  1-foot  barrier,  but  the  solutions 
for  the  8-foot  case  are  very  close  down  to  kR  •  10. 

Comparing  Figure  4.21  with  Figure  4.19  shows  the  expected  result 
that  the  higher  barrier  provides  the  greater  attenuation. 
Mathematically,  this  is  a  result  of  the  diffraction  coefficient  taking 

t 

smaller  values  as  the  receiver  moves  deeper  into  the  shadow  zone.  For 
soft  grounds,  however,  a  counteracting  effect  takes  place  as  the 
barrier  height  increases.  The  reflection  angles  of  the  source-ground- 
barrier  path  and  the  barrier-ground-receiver  path  Increase,  and  hence 
the  "strength"  of  the  ground-reflected  rays  gets  larger.  The  net 
result  can  be  seen  most  clearly  when  Figure  4.18  is  compared  to  Figure 
4.20  (attenuation  or  insertion  loss).  The  attenuation  is  higher  for 
the  8-foot  barrier  for  kR  less  than  about  50,  but  the  attenuation  is 
generally  lower  for  kR  above  that.  Apparently,  then,  the  "diffraction 
effect"  is  the  primary  factor  at  the  lower  values  of  kR,  while  the 
"Impedance  effect"  dominates  at  the  higher  values. 

When  both  the  source  and  receiver  are  raised  above  the  ground  to 
heights  of  7.2  feet  and  7.1  feet,  respectively,  the  data  plotted  in 
Figure  4.22  result.  Since  the  geometry  chosen  here  defines 
diffraction  angles  and  <p  equal  to  those  for  the  case  in  Figure 
4.18,  the  differences  between  the  two  sets  of  data  are  due  entirely  to 
the  presence  of  the  ground  (the  half-plane  curve  is  identical  for  the 


!  Calculated  Data  for  Attenuation  as  a  Function  of  kR.  Q-term,  Plane  Wave,  and 
Solutions.  Large  Source  and  Receiver  Heights,  Large  Barrier  Height,  Low  Impedanc 
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two  geometries).  The  barrier  Is  at  8.0  feet  and  the  ground  Impedance 
Is  Z.,  -  (1.5,  2.0).  The  first  i .  k  in  the  attenuation  is  due  to  the 
ground  impedance;  it  occurs  at  a  low  value  of  kR  because  the 
reflection  angles  are  now  large  (recall  the  trend  noted  in  Section 
4.1.1  that  as  iff  increased,  the  peak  shifted  to  lower  kR).  The  series 
of  peaks  and  troughs  for  kR  _>  10  are  due  to  path  length  interferences 
among  the  four  ray  paths  considered  in  the  Edge-Plus-Images  model. 
Besides  the  fact  that  the  plane  wave  solution  is  valid  here  for  kR 
larger  than  10,  the  interesting  feature  of  the  graph  is  that  the  half¬ 
plane  solution  seems  to  predict  the  "average"  attenuation. 

The  next  series  of  graphs,  Figures  4.23  -  4.24  show  the  effects 
of  having  different  impedances  on  each  side  of  the  barrier.  For  all 
four  cases,  the  source  height  is  0.25  feet;  the  receiver  height  is  8.0 
feet;  the  source  and  receiver  are  60.0  feet  from  the  barrier;  and  the 
barrier  is  15.0  feet  high.  The  geometry  is  realistic  for  a  practical 
highway  noise  barrier  problem.  First,  the  data  for  a  hard  (ZN  » 
[100.0,  150.0])  ground  beneath  both  source  and  receiver  is  plotted  in 
Figure  4.23  (a)  (plane  wave  and  Q-term  solutions  only).  For  such  a 
high  impedance,  both  solutions  are  accurate.  When  the  ground  on  the 
receiver  side  is  assigned  a  constant  impedance  Z£  *  (2.0,  2.5),  the 
data  in  Figure  4.23  (b)  result.  The  overall  attenuation  is  increased 
at  lower  frequencies  as  the  first  peak  in  Figure  4.23  (a)  shifts  to 
the  left;  also,  the  solutions  begin  to  show  differences.  The  effect 
of  the  impedance  at  high  frequencies  is  simply  to  "smooth"  the  peaks 


an  troughs,  as  the  average  attenuation  changes  little. 
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If  Che  soft  ground  were  located  under  the  source  instead  of  the 
receiver,  Figure  4.24  (a)  would  describe  the  resulting  attenuation. 
This  shows  several  noteworthy  features.  First,  the  entire  curve  is 
shifted  up  relative  to  the  hard-ground  case  in  Figure  4.23  (a). 

Second,  the  first  peak  of  attenuation  centered  at  about  kR-100  does 
not  shift;  and  third,  there  is  no  smoothing  of  the  high-frequency 
peaks  (although  the  average  is  now  higher).  All  of  these  results  are 
due  to  the  fact  that,  unlike  the  receiver,  the  source  is  very  close 
to  the  ground  and  the  propagation  is  at  near-grazing  angles,  resulting 
in  higher  losses.  Since  the  ’’incident"  field  at  the  barrier  edge  is 
reduced  by  the  source  ground  impedance  layer,  the  entire  pattern  of 
receiver-side  Interference  is  shifted  upward  but  not  changed  in 
character.  Finally,  when  the  ground  on  both  sides  is  soft,  as  in 
Figure  4.24(b),  the  effects  are  cumulative:  an  overall  increase  in 
attenuation,  a  shifting  of  the  first  peak,  and  a  smoothing  of  the 
Interference  peaks. 

To  conclude  this  section,  a  case  for  which  the  incident  ray 
strikes  the  barrier  edge  obliquely  will  be  investigated.  The  geometry 
is  the  same  as  in  the  previous  case  except  that  Zq  >100  (Figure  4.17 
[b]),  giving  an  angle  0  ■  50°.  The  ground  impedance  is  ZN  "  <100-0’ 
150.0).  Figure  4.25  shows  the  resulting  attenuation,  with  the  results 
of  Figure  4.23  (a)  superimposed  for  comparison.  Since  a  factor  of 
sin0  appears  in  the  denominator  of  the  diffraction  coefficient  given 
in  Equation  (3.4),  the  coefficient  itself  increases  as  the  angle  3 
decreases.  The  numerical  predictions  thus  indicate  a  lower 


mparison  of  Results  for  Perpendicular  and  Oblique  Incidence  on  the  Edge  as 
Function  of  kR.  Q-term,  Plane  Wave,  and  Half-Plane  Solutions.  Hard 
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attenuation  for  the  oblique  Incidence  case.  As  this  Is  not  a  ground- 
related  effect,  the  half-plane  solution  (long-dashed  line)  predicts 
the  same  trend. 

The  geometry  and  impedance  conditions  for  the  plots  presented  in 
this  section  are  summarized  in  Table  4.6. 

4.2.2  The  Dependence  on  Barrier  Height 

Half-plane  diffraction  theory  predicts  the  simple  formula  for 
fixed  source  and  receiver  positions:  the  greater  the  barrier  height, 
the  greater  the  attenuation.  While  this  is  generally  true  in 
practice,  consideration  of  the  ground  reflections  and  impedance 
conditions  greatly  affects  the  predicted  attenuation,  and  sometimes 
leads  to  unexpected  results. 

The  data  plotted  in  Figures  4.26  and  4.27  compare  the  predicted 
attenuation  and  insertion  loss  for  a  soft  (ZN«[1.5,  2.0))  ground  with 
that  for  a  relatively  hard  (Z^*[15.0,  20.0])  ground.  The  source  and 
reciver  are  close  to  the  surface  (0.5  feet)  and  relatively  far  from 
the  barrier  (60  feet),  and  kR  is  assigned  a  small  value  (5.0).  The 
barrier  height  ranges  from  0.1  to  25  feet.  Several  observations  can 
be  made.  First,  the  low  value  of  kR  is  responsible  for  the  large 
differences  between  the  plane  wave  solution  and  the  Q-term  solution. 
Even  when  the  barrier  is  at  maximum  height  (the  ground  reflection 
angles  are  near  25°),  there  remains  a  7  dB  difference  for  the  low- 
impedance  case  and  a  1  dfe  difference  for  the  higher-impedance  case. 
Second,  the  overall  effects  of  diffraction  are  minimal  on  the  Q-term 
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predictions  as  the  insertion  loss  is  close  to  zero  over  most  of  the 
range  of  barrier  height.  This  is  simply  a  consequence  of  the  very  low 
frequencies  being  considered  here. 

Figures  4.28  and  4.29  show  data  for  the  same  geometry  as  above, 
but  here  the  value  of  kR  has  been  increased  by  a  factor  of  ten  to 
kR*50.  The  plane  wave  and  Q-term  solutions  are  more  in  accord  here. 

In  fact,  for  barrier  heights  above  12  feet  the  predictions  are  within 
1  dB  for  either  impedance  condition.  The  fact  that  a  broad  peak  in 
attenuation  appears  for  the  soft  ground  and  not  for  harder  ground 
indicates  that  this  is  an  “impedance  effect"  similar  to  that  noticed 
in  previous  ground  propagation  curves. 

The  two  insertion  loss  curves  in  Figures  4.28  (b)  and  4.29  (b) 
reveal  an  interesting  trend  that  was  noted  in  the  previous  section. 

That  is,  the  insertion  loss  is  near  zero  for  high  barriers  on  the  soft 
ground,  while  it  steadily  Increases  with  height  for  barriers  on  the 
harder  ground.  In  fact,  for  the  soft  ground,  the  insertion  loss  shows 
a  steady  decline  (Figure  4.28  [b])  for  heights  above  4  feet; 
therefore,  increasing  the  barrier  height  beyond  this  value  would  be 
counterproductive  (in  the  noise  control  sense). 

The  counteracting  Impedance  and  diffraction  effects  mentioned  in 
the  previous  section,  are  apparent  in  the  next  pair  of  graphs.  In 
Figure  4.30  (a)  and  (b),  the  source  and  receiver  have  been  raised 
above  the  ground  to  heights  of  4  feet,  but  all  other  conditions  are 
the  same  as  in  the  previous  curves.  For  the  case  where  »  (1.5, 

2.0),  the  attenuation  stays  relatively  constant  as  the  barrier  height 
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Figure  4.29  Calculated  Data  for  (a)  Attenuation  and  (b)  Insertion 
Loss  as  a  Function  of  Barrier  Height.  Q-term,  Plane 
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Increases.  On  the  other  hand,  each  ray  path  in  the  Edge-Plus-Images 
model  (Figure  3.4)  should  experience  a  trend  similar  to  that  shown  by 
the  single-ray  half-plane  diffraction  solution  plotted  on  the  same 
graph.  That  is,  each  of  the  four  paths  is  attenuated  more  heavily  as 
the  barrier  height  increases.  However,  the  strength  of  each  ground- 
reflected  path  gets  larger  as  the  reflection  angle  Increases 
(approximately,  from  4°  to  25°  over  the  range  of  barrier  heights). 

This  effect  offsets  the  Increase  in  attenuation  due  to  diffraction. 
Finally,  it  can  be  seen  that  this  phenomenon  does  not  occur  for  the 
harder  (ZN-[100.0,  150.0])  ground  shown  in  Figure  4.30  (b);  the 
impedance  effect  is  not  pronounced  since  the  ground  appears  hard 
throughout  the  range  of  reflection  angles.  Thus,  the  attenuation 
steadily  increases  with  barrier  height. 

To  conclude  this  section,  two  additional  sets  of  curves  are 
presented  in  Figures  4.31  and  4.32  for  large  values  of  kR.  The 
parameter  values  are  indicated  on  the  graphs.  Essentially,  both  plots 
are  for  soft  ground  but  the  former  has  the  receiver  well  above  the 
ground  while  the  latter  has  it  very  close  to  the  ground.  The  plane 
wave  solution  is  accurate  in  either  case  over  the  full  range  of 
barrier  heights.  Other  important  observations  are  that  path  length 
interference  effects  are  present  for  the  elevated  receiver,  and  that 
the  insertion  loss  is  negative  when  the  source  and  receiver  are  near 
the  ground. 

Table  4.7  presents  a  summary  of  the  geometrical  and  impedance 
conditions  for  the  numerical  results  plotted  in  this  section. 
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Figure  4.31  Calculated  Data  for  (a)  Attenuation  and  (b)  Insertion 
Loss  as  a  Function  of  Barrier  Height.  Q-term,  Plane 
Wave,  and  Half-Plane  Solutions.  Sht*0.1,  Rht-8.0, 
kR-1000,  and  Z  -  (1.3,  1.2). 


Figure  4.32  Calculated  Data  for  (a)  Attenuation  and  (b)  Insertion 
Loss  as  a  Function  of  Barrier  Height.  Q-term,  Plane 
Wave,  and  Half-Plane  Solutions.  Sht*0.5,  Rht*0.5, 
kR*400,  and  ZM  ■  (1.5,  2.0). 


Table  4.7 


4.2.3  Dependence  on  Diffraction  Angle  4> 


The  diffraction  angle  $  defined  in  Figure  4.17  naturally  depends 
on  the  barrier  height  and  the  receiver  position.  Consequently,  many 
observations  regarding  the  dependence  of  the  field  on  4>  have  already 
been,  or  can  be,  made  from  the  numerical  data  presented  thus  far. 
However,  for  completeness,  several  cases  will  be  presented  here 
showing  this  dependence  explicitly. 

The  numerical  data  predicted  from  the  plane  wave,  Q-term,  and 
half -plane * theories  for  a  source  0.2  feet  above  the  ground  and  a 
barrier  of  height  5.0  feet  are  plotted  in  Figures  4.33  -  4.36. 
Sequentially,  the  graphs  show  a  small  kR  coupled  with  a  low  or  a  high 
ground  Impedance,  and  a  large  kR  with  a  low  or  a  high  ground 
impedance.  In  these  plots,  the  receiver  distance  from  the  edge,  R^, 
is  held  constant  as  the  angle  <J>  is  varied.  Consequently,  the  receiver 
revolves  about  the  the  edge  point  E,  moving  from  its  starting  point  on 
the  ground  behind  the  barrier  (<f>  maximum)  to  its  final  position 
directly  above  the  edge  ($  ■  0°).  At  some  point  in  this  angular 
range,  the  receiver  will  cross  the  incident  shadow  boundary  and  pass 
into  the  illuminated  region  of  the  source.  For  the  geometry 
considered  here,  the  shadow  boundary  from  the  source  occurs  at  69.2°, 
and  this  value  (^g^)  *8  narked  by  a  vertical  line  on  the  graphs. 


The  values  marked  on  the  axis  for  the  angle  f  are  actually 
absolute  values;  the  cylindrical  coordinate  system  used  in  the 
development  of  the  model  ieflne  the  angle  $  to  be  negative  on  the 
receiver  side  of  the  barrier. 
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Figure  4.33  Calculated  Data  for  (a)  Attenuation  and  (b)  Insertion 

Loss  as  a  Function  of  Diffraction  Angle.  Q-term,  Plane 
Wave,  and  Half-Plane  Solutions  Sht-0.2,  *12.0, 

kR* 5.0,  Bht-5.0,  and  Z  -(1.5,  2.0). 
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Figure  4.34  Calculated  Data  for  (a)  Attenuation  and  (b)  Insertion 

Loss  as  a  Function  of  Diffraction  Angle.  Q-term,  Plane 
Wave,  and  Half-Plane  Solutions.  Sht=0.2,  R_=12.0, 
kR»5.0,  Bht”5.0,  and  Z  =(15.0,  20.0). 
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Figure  4.35  Calculated  Data  for  (a)  Attenuation  and  (b> 
Loss  as  a  Function  of  Diffraction  Angle.  Q-terra. 
Wave,  and  Half-Plane  Solutions.  Sht-0.2,  R  ,  * ;  2  . 
Bht-5.0,  and  Zvt-(1.5,  2.0). 
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Figure  4.36  Calculated  Data  for  (a)  Attenuation  and  (b)  Insertion 
Loss  as  a  Function  of  Diffraction  Angle.  Q-term,  Plane 
Wave,  and  Half-Plane  Solutions.  Sht-0.2,  R.-12.0,  kit-100, 
Bht-5.0,  and  Z^-(15.0,  20.0). 
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The  only  significant  differences  between  the  plane  wave  and  the 
Q-term  predictions  (for  this  particular  barrier  height)  are  for  the 
low-impedance  (Z^-[1.5,2.0]),  small  kR  (kR*5.0)  case  shown  in  Figure 
4.33.  For  the  distances  considered  here  (R-24.0),  this  value  of  kR 
corresponds  to  a  frequency  of  37  Hz.,  and  it  is  very  unlikely  that  any 
practical  material  (in  air)  would  possess  such  a  low  impedance  at  this 
low  frequency.  The  cases  with  a  low  value  of  kR  show  a  broad  peak  in 
the  attenuation  or  insertion  loss  response,  whereas  the  cases  with  the 
higher  value  of  kR  show  the  expected  steady  decline  as  the  receiver 
approaches  the  incident  shadow  boundary.  It  is  interesting  .to  note 
that  for  the  large  kR  case,  the  simple  half-plane  diffraction  model 
overpredicts  the  attenuation  for  the  hard  ground  and  underpredicts  it 
for  the  soft  ground. 

Some  further  comments  about  the  shadow  boundary  may  be 
appropriate  here.  It  is  well  known  that  the  presence  of  the  Fresnel 
integrals  in  the  diffraction  coefficient  (Equation  13.4]),  ensures 
that  the  diffracted  field  is  continuous  across  the  shadow  boundary. 

The  numerical  data  plotted  here  confirm  this  fact;  no  abrupt  changes 
occur  in  either  the  attenuation  or  insertion  loss  data  in  the  vicinity 
of  the  shadow  boundary.  It  is  also  well  known  that  the  diffracted 
field  at  the  shadow  boundary  itself  (for  large  kR)  assumes  a  value 
one-half  that  of  the  incident  field  there.  This,  of  course,  applies 
to  the  diffraction  of  a  single  ray  only.  Therefore,  the  data  in 
Figures  4.35  (a)  and  4.36  (a)  show  a  6  dB  attenuation  for  the  half¬ 
plane  predictions  but  a  different  value  when  the  ground  is  taken  into 
account  (the  plane-wave  or  Q-term  data). 
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As  mentioned  In  the  description  of  the  Edge-Plus-Images  model  in 
Chapter  3,  the  source  image  in  the  ground  also  defines  its  own  shadow 


boundary*  That  is,  there  is  a  "shadow”  and  an  "illuminated”  region 
for  the  reflected  acoustic  energy  as  well  as  for  the  direct  energy 
from  the  source.  As  the  angle  $  decreases,  the  receiver  may  also 
cross  this  "image  shadow  boundary,"  but,  again,  there  should  be  no 
discontinuity  in  the  predicted  levels.  Thiti.  xs  again  confirmed  by  the 
data,  since  the  shadow  boundary  for  the  image  source  occurs  within  a 
few  degrees  of  the  Incident  shadow  boundary  shown  on  the  graphs. 

The  relevant  parameters  for  the  plots  in  this  section  are 
summarized  in  Table  4.8. 

4.2.4  Received  Spectra  for  a  Practical  Ground  Cover 

As  noted  in  Section  4.1.3,  the  ground  impedance  is  rarely 
constant  over  a  wide  range  of  frequency  values.  Therefore,  as  a 
practical  example,  the  Edge-Plus- Images  prediction  model  has  been 
applied  to  a  barrier  on  a  ground  plane  characterized  by  the  third- 
octave  values  of  Impedance  given  in  Table  4.3.  The  source  and 
receiver  are  very  close  to  the  ground  (0.25  feet),  and  the  barrier 
height  is  15.0  feet. 

The  data  plotted  in  Figure  4.37  are  for  a  perfectly  rigid  ground. 
Of  course,  the  plane  wave  and  Q-term  solutions  are  Identical,  both 
indicating  an  attenuation  that  gradually  increases  with  frequency 
(diffraction  coefficient  decreasing  with  increasing  kR)  until  a  sharp 
peak  due  to  ground  reflection  Interference  (path  length  difference)  is 
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Figure  4.37  Calcualted  Data  for  (a)  Attenuation  and  (b)  Insertion 
Loss  as  a  Function  of  Frequency.  Q-tenn  and  Plane  Wave 
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reached.  The  Insertion  loss  follows  this  trend  exactly.  When  the 
ground  is  grass-covered,  the  data  in  Figure  4.38  is  generated.  An 
important  feature  here  is  that  the  plane  wave  solution  remains  valid 
even  for  very  low  kR.  This  can  be  explained  with  reference  to  Table 
4.3,  where  it  can  be  seen  that  for  low  frequencies,  the  impedance  is 
rather  high.  This  tends  to  counteract  any  "low  kR"  differences 
between  the  plane  wave  and  Q-term  solutions. 

As  has  been  observed  in  previous  cases,  the  attenuation  shows  a 
broad  peak  in  the  response  for  the  impedance-covered  ground.  However, 
Figure  4.38  (b)  shows  that  the  corresponding  insertion  loss  drops 
rapidly  in  this  region.  Again,  the  barrier  is  "preventing"  the  large 
grazing-incidence  propagation  loss  to  take  place  over  the  soft  ground. 

As  a  final  exercise,  a  case  has  been  analyzed  in  which 
significant  differences  do  exist  between  the  plane  wave  and  Q-term 
predictions,  in  spite  of  the  high  impedances  at  low  kR  values.  That 
is,  the  barrier  height  has  been  reduced  (Bht-0.5  feet)  to  force  near- 
grazing  incidence  for  the  ground  reflected  rays.  The  numerical  data 
are  plotted  in  Figure  4.39  (a)  and  (b),  where  the  plane  wave  solution 
is  seen  to  be  up  to  30  dB  Inaccurate.  This  configuration  of  source, 
receiver,  and  barrier  is  certainly  not  representative  of  noise  control 
applications,  but  it  is  Interesting  to  note  the  single  sharp  peak  in 
the  Insertion  loss  spectra. 

Data  for  a  variety  of  cases  have  been  presented  in  this  chapter, 
mainly  to  illustrate  some  key  differences  between  the  classical  plane 
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Figure  4.39  Calcualted  Data  for  (a)  Attenuation  and  (b)  Insertion 
Loss  as  a  Function  of  Frequency.  Q-term  and  Plane  Wave 
Solutions.  Sht«0.25,  Rht-0.25,  Bht«0.5,  and  Z^-GRASS. 


wave  solution  to  the  propagation  problem  and  the  more  exact  spherical 

wave  solution.  A  secondary  motivation  was  to  collect  a  large  sample 

* 

of  data  for  future  reference.  In  view  of  the  many  parameters 
involved,  this  same  data  could  provide  much  more  information  than  has 
been  presented  here.  Furthermore,  numerical  results  for  either  the 
ground  propagation  solution  or  the  Edge-Plus-Images  model  are 
relatively  easy  to  generate  with  a  digital  computer,  and  so  a 
particular  feature  or  trend  in  any  one  plot  could  be  Investigated  in 
detail,  if  desired. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

The  primary  objective  of  the  research  behind  this  study  was  to 
derive  a  more  exact  analytical  solution  to  the  problem  of  point  source 
propagation  over  an  impedance  plane.  This  goal  has  been  met  in 
obtaining  the  asymptotic  solution  described  in  Chapter  2.  A  second 
objective  was  to  incorporate  the  new  solution  into  a  barrier  model 
that  would  account  for  ground  reflections  in  addition  to  diffraction. 
This  was  also  accomplished,  and  the  so-called  Edge-Plus- Images  model 
discussed  in  Chapter  3  was  developed. 


The  mathematical  problem  was  formulated  as  a  boundary  value 
problem,  in  terms  of  the  acoustic  Helmholtz  equation  in  cylindrical 
coordinates  (Equation  [2.18])  and  the  local  reaction  boundary 
condltiolis  (Equation  [2.40]).  Using  Hankel  Transforms,  an  integral 
solution  was  obtained  (Equation  [2.48])  which  contained  a  pole  and  a 
branch  point.  The  pole  was  replaced  by  an  equivalent  Integral 
(Equation  [2.51]),  and  the  resulting  double  Integral  was  found  to 
contain  a  term  resembling  the  Sommerfeld  representation  for  a  point 
source.  The  double  integral  thus  reduced  to  the  single  Integral  given 
by  Equation  (2.57). 

The  major  part  of  the  research  effort  was  directed  toward 
evaluating  the  integral  in  Equation  (2.57).  Other  researchers  have 
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arrived  at  similar  Integrals  when  studying  this  problem  an)  have 
invoked  such  asymptotic  techniques  as  the  method  of  steepest  descents 
for  its  evaluation.  The  approach  used  in  this  study  was  to  apply  the 
variable  transformation  given  by  Equation  (2.60)  to  the  integrand.  In 
this  manner,  the  complicated  exponential  was  reduced  to  a  relatively 
simple  polynomial  expression.  The  remainder  of  the  integrand  was 
transformed  using  standard  techniques,  and  the  final  real-axis 
integral  given  in  Equation  (2.93)  was  obtained. 

Since  the  transformed  integral  is  still  intractable  due  to  the 
square  root  function  in  the  denominator,  a  Taylor  series  expansion  of 
the  latter  was  sought  and  obtained  (Equation  [2.95]).  The  important 
result  here  is  that  a  general  expression  for  the  nth  Taylor 
coefficient  was  derived,  making  available  (computatlonallv)  as  many 
terms  in  the  series  as  might  be  desired  or  necessary.  This  expansion 
of  the  denominator  allowed  the  original  integral  in  Equation  (2.93)  to 
be  written  as  an  infinite  sum  of  integrals  as  in  Equation  (2.105). 

Tb'  sum  is  an  asymptotic  series  in  the  parameter  "kR^.”  , 

The  final  step  in  the  formal  derivation  was  to  perform  the  term- 
by-term  integration,  that  is,  to  evaluate  the  component  integrals  in 
the  asymptotic  series.  Each  integral  was  resolved  in  terms  of 
parabolic  cylinder  functions  (Equation  [2.106]),  which  in  turn  were 
expressed  in  terms  of  the  complex  complementary  error  function.  From 
the  perspective  of  this  study,  the  most  important  feature  of  the 
parabolic  cylinder  functions  was  found  to  be  their  recursive  nature. 
Specifically,  the  recursion  relations  for  the  parabolic  cylinder 
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functions  led  directly  to  recursion  relations  for  the  integral  terms 
in  the  asymptotic  series  (Equations  [2.117]  -  [2.119]).  Thus,  the 
final  result  of  the  mathematical  derivation,  the  closed-form 
asymptotic  series  given  by  Equation  (2.123)  or  (2.128),  was  obtained. 

The  formal  solution  was  then  re-phrased  in  several  useful  forms. 
The  "Q-term”  form  in  Equation  (2.129),  with  Q  defined  by  Equation 
(2.135),  lends  a  physical  Interpretation.  Namely,  the  quantity  Q 
could  be  thought  of  as  a  "spherical  wave  reflection  coefficient," 
since  the  equation  in  this  form  is  analogous  to  that  employing  the 
classical  plane  wave  reflection  coefficient  (Equation  [2.5]).  The  Q- 
term  form  of  solution  is  preferred  here;  for  reference,  it  is  re¬ 
stated  in  Figure  5.1.  The  ”F-term”  form  given  by  Equation  (2.137), 
with  F  defined  by  Equation  (2.141),  presents  the  solution  as  one 
containing  a  "correction"  term  for  the  classical  plane  wave  solution. 
The  form  of  the  solution  given  in  Section  2. 3. 5. 4  shows  that  the 
Hankel  function  can  be  extracted  from  the  asymptotic  series  (other 
researchers  have  generated  approximate  solutions  containing  the  Hankel 
function).  Finally,  a  form  of  solution  using  only  the  first  term  in 
the  asymptotic  series  was  presented  in  Section  2. 3. 5. 5. 

For  the  case  of  perpendicular  incidence  (sin<|>  ■  1.0),  an  exact 
solution  was  derived  in  terms  of  the  exponential  integral  (Equation 
[2.158]).  Small-  and  large-argument  expansions  for  the  exponential 
integral  were  available  and  are  given  by  Equations  (2.156)  and 
(2.157),  respectively.  It  was  found  that  the  solution  using  the  first 
term  in  the  latter  expansion  is  identical  to  the  plane  wave  solution 


for  this  case 
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The  next  phase  of  the  research  was  directed  toward  incorporating 
the  ground  propagation  solution  into  a  barrier  diffraction  model. 
Following  a  brief  review  of  the  underlying  diffraction  concepts  in 
Section  3.1,  the  principal  result  of  this  phase — the  so-called  Edge- 
Plus-Images  model — is  described  in  Section  3.2.  The  model  is 
constructed  using  four  separate  half-plane  diffracted  ray  paths.  That 
is,  a  direct-diffracted  ray,  a  reflected-diffracted  ray,  a  diffracted- 
reflected  ray,  and  a  reflected-diffracted-reflected  ray  (Figure  3.4) 
are  combined  coherently,  and  the  ground-reflected  rays  are 
appropriately  modified  by  the  spherical  wave  reflection  coefficient. 

A  large  quantity  of  numerical  data  generated  from  the  ground 

l 

propagation  solution  and  from  the  Edge-Plus-Images  diffraction  model 
was  analyzed  next.  The  dependence  of  the  predictions  on  several 
important  parameters  was  investigated  in  Chapter  4  in  a  comprehensive 
series  of  graphs.  For  the  ground  propagation  section,  the  dependence 
on  the  parameter  kR,  the  reflection  angle  ly,  and  the  ground  impedance 
ZN  was  studied,  and  for  the  barrier  section,  the  influence  of  barrier 
height  and  diffraction  angle  $  was  studied  in  addition. 

5.2  Conclusions 

In  the  discussion  of  the  numerical  results  from  both  the 
propagation  solution  and  the  barrier  diffraction  model,  many  specific 
observations  of  trends  and  features  in  the  data  were  pointed  out. 

Here,  several  general  conclusions  will  be  drawn  and  discussed. 
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The  asymptotic  series  solution  is  very  accurate.  The  only 
noticeable  difference  between  the  data  predicted  from  the  asymptotic 
series  solution  and  the  data  from  the  numerical  integration  of  the 
exact  solution  occurred  at  very  low  values  of  kR.  Deviations  in  this 
region  are  expected,  but  the  surprising  fact  is  that  these  deviations 
cease  to  exist  for  values  of  kR  that  are  still  relatively  small.  For 
example.  Figure  4.1  shows  no  difference  between  the  exact  and 
asymptotic  solutions  for  kR  down  to  a  value  kR-3.  And  Figure  4.2 
shows  only  0.25  difference  for  kR*0.1.  Re-examining  the  integrand 
trends  presented  in  Figure  2.7  (a),  (b),  and  (c)  might  provide  some 
insight  into  why  this  is  possible.  For  the  low  kR  values  cited  here, 
the  integrand  will  show  significant  amplitudes  outside  the  radius  of 
curvature  of  the  Taylor  series  (t  ).  Presumably,  then,  the 
oscillatory  nature  of  the  integrand  prevents  any  net  contribution  to 
the  integral  for  all  but  the  smallest  values  of  t  (for  which  the 
Taylor  series  is  convergent). 

The  first  term  in  the  asymptotic  series  usually  provides 
sufficient  accuracy.  For  the  geometries  and  impedance  conditions 
considered  in  this  study,  the  maximum  deviation  between  the  F-term 
(that  is,  from  Equation  [2.150])  solution  and  the  full  series 
solution  can  be  seen  in  Figure  4.3.  Here,  a  3  dB  deviation  exists 
over  a  wide  range  of  kR.  However,  this  is  an  extreme  case  (grazing 
incidence,  very  soft  ground),  and  a  more  typical  behavior  (Figures  4.5 
and  4.6,  for  example)  shows  agreement  with  the  exact  solution  for 


larger  kR  and  only  one  or  two  dB  deviation  for  small  kR.  It  is 
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interesting  to  note  that  for  extremely  small  kR,  the  solution  using 
only  the  first  term  in  the  asymptotic  series  may  give  more  accuracy 
than  the  solution  using  more  terms— these  terms  are  probably  diverging 
here— as  can  be  seen  in  Figures  4.1  (b)  or  4.3,  for  example. 

The  magnitude  of  the  difference  between  the  solution  using  the 
first  term  only  and  that  using  the  full  series  is  highly  dependent  on 
the  .values  of  the  parameters  kR,  Z^,  and  .  The  effects  of  these 
parameters  are  inter-related,  and  it  is  difficult  to  determine  which 
is  the  “controlling"  factor.  The  ground  impedance  seems  to  have  a 
very  strong  Influence.  That  is,  even  for  relatively  large  kR  and 
large  ijs,  a  small  value  for  Z  will  still  produce  deviations  (for 
instance,  in  Figure  4.8  [b]). 

The  simple  plane  wave  solution  may  be  accurate  over  a  wide  range 
of  parameter  values.  It  is  definitely  Inaccurate  for  small  values  of 
all  three  parameters  and  definitely  accurate  for  large  values  of  all 
three.  No  general  rule  can  be  stated  for  intermediate  values,  yet  the 
numerical  data  in  Chapter  4  may  serve  as  a  guide. 

The  ground  reflections  must  be  taken  into  account  when 
considering  diffraction  by  a  barrier.  For  all  the  cases  considered 
here,  the  half-plane  prediction  curve  (diffraction  alone)  was 
significantly  inaccurate,  even  for  very  practical  geometries  and 
impedance  conditions.  The  only  possible  exception  to  this  conclusion 
is  the  case  where  the  source  and  receiver  are  well  above  the  ground, 
for  which  the  half-plane  theory  seems  to  predict  the  "average” 
attenuation  through  the  interference  maxima  and  minima  (Figure  4. 25). 


barrier  la  present.  This  la  bacauaa  of  the  larger  ground  relfectlon 
angles.  The  Edge-Plue-Inagea  aodel  ssy  only  require  the  uae  of  the 
more-exact  spherical  wave  reflection  coefficient  for  very  short 
barrier  heights,  source  and  receiver  near  the  ground,  and  a  relatively 
low  ground  impedance.  Practical  cases  (in  outdoor  noise  control,  at 
least)  show  very  little  difference  between  the  solutions  (Figures  4.43 
and  4.44). 


testlons  for  Further  Research 


From  the  mathematical  standpoint,  more  rigor  could  be  applied  to 
several  points  in  the  derivation  of  the  ground  propagation  solution. 

An  investigation  of  the  convergence  of  the  asymptotic  series,  and  a 
more  thorough  treatment  of  the  branch  cuts  and  their  dependence  on  the 
angle  and  the  admittance  could  be  made.  For  the  Edge-Plus-Images 
model,  the  diffraction  coefficient  could  be  adapted  to  other  types  of 
barriers,  such  as  wedges,  trapezoids,  or  parallel  edges. 

An  expanded  series  of  sensitivity  tests  could  also  be  undertaken 
to  determine  more  precise  relationships  among  the  parameters.  In 
addition,  numerical  data  for  the  absorptive  barrier  might  be  generated 
and  compared  to  the  rigid  barrier. 

Finally,  it  may  be  interesting  to  apply  the  theoretical  results 
of  the  present  study  to  other  fields  such  as  underwater  acoustics, 
seismology,  or  electromagnetics,  for  which  certain  values  of  the 
parameters  kR  and  ZN  considered  "unrealistic”  in  air  acoustics  might 


be  more  natural 
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